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COMPUTABLE  REAL  ANALYSIS 


B.  J.  MacLennan 

Abstract: 

We  present  a  model  of  the  real  numbers  that  is  completely  finitutic.  Every  real  number  is  represented 
by  a  finite  structure  —  specifically,  a  finite  sequence  of  symbols  from  a  finite  alphabet.  All  of  the  arith- 
metic  operations  on  reals  are  also  finite  and  can  be  evaluated  on  a  computer.  We  allow  nothing  that  can¬ 
not  >e  described  by  a  finite  algorithm  -  whether  numbers  or  sequences  or  functions.  This  development  is 
carried  through  the  fundamental  theorem  of  calculus. 

1.  Introduction 

In  this  report  we  develop  systematically  the  foundations  of  a  theory  of  computable  real  analysts.  This 
is  based  on  the  view  that  the  only  objects  that  it  makes  sense  to  talk  about  are  finite  objects,  and  that  the 
only  operations  it  makes  sense  to  talk  about  are  those  that  can  be  completed  in  a  finite  number  of  opera¬ 
tions. 

This  approach  does  not  exclude  all  use  of  notions  of  infinity.  In  particular,  it  permits  potential  infini¬ 
te*,  provided  the  operation  generating  the  infinity  can  be  finitely  described,  or  is  itself  the  result  of  a  pro¬ 
cess  that  can  be  finitely  described,  and  so  forth,  so  that  it  can  be  reduced  to  a  finite  process  in  a  finite 
number  of  steps.  For  example,  we  will  permit  the  use  of  infinite  sequences,  provided  that  their  enumera¬ 
tion  functions  are  finitely  describable  and  finitely  computable.  We  eschew  all  use  of  actual  infinitie*,  that 
is,  all  infinities  that  are  irreducible  to  finite  objects. 

This  approach  is  most  easily  accomplished  by  restricting  our  attention  to  finite  objects  and  computable 
functions.  In  other  words,  we  will  restrict  our  attention  to  computer  programs  operating  on  computer 
data.  Computer  programs  and  data  are  both  finite  objects,  but  they  provide  a  convenient  means  for  deal¬ 
ing  with  potential  infinities. 

Notice  that  in  this  approach  every  real  number  is  represented  by  a  finite  structure  -  that  is,  a  finite 
string  of  symbols  chosen  from  a  finite  alphabet.  In  particular,  irrational  numbers  are  represented  finitely. 
This  reduction  to  finite  structures  is  accomplished  by  identifying  real  numbers  with  programs  that  com- 
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pute  sequences  of  rational  approximations  to  the  intended  real  number. 

Various  operations  and  functions  of  real  numbers  (e.g.  the  extracting  of  limits  of  sequences  of  reals)  are 
also  defined  as  programs,  which  means  that  these  operations  and  functions  are  computable  in  a  finite 
number  of  operations.  Thus  we  have  the  foundations  of  a  completely  finitiitic  theory  of  real  analysis. 

There  are  many  models  of  computation  that  could  be  used  in  this  development.  There  are  several  rea¬ 
sons  why  we  have  chosen  to  use  functional  programming  as  a  basis  for  our  development  of  the  foundations 
of  real  analysis: 

1.  Functional  programming  is  mathematical  in  its  style.  Thus  it  lends  itself  to  the  definition  of 
mathematical  ideas  and  is  conducive  to  proofs  of  mathematical  properties. 

2.  Functional  programming  provides  high-level,  powerful  means  for  manipulating  computable  functions. 
This  is  especially  helpful  in  defining  operations  on  potentially  infinite  objects. 

3.  Functional  programs  can  be  straight-forwardly  reduced  to  formulas  in  the  lambda  calculus,  a  well- 
known  model  of  computation.  This  reduction  founds  real  analysis  on  a  small  number  of  basic  ideas. 

A  secondary  goal  of  this  report  is  to  demonstrate  the  use  of  functional  programming  concepts  in  construc¬ 
tivist  mathematics. 

We  have  omitted  many  of  the  proofs  that  our  definitions  are  correct,  or  that  they  satisfy  expected  pro¬ 
perties  (commutativity,  associativity,  etc.).  In  most  cases  these  proofs  are  routine  and  will  not  be  missed. 
However,  we  have  frequently  given  informal  demonstrations  that  our  definitions  are  correct. 

In  this  report  we  follow  closely  the  development  in  Errett  Bishop’s  Foundation i  of  Conttructivc 
Analgiii  (McGraw-Hill,  New  York,  1967).  Both  the  definitions  and  order  of  presentation  have  been  based 
on  his,  although  they  have  been  adapted  where  necessary  to  the  functional  approach.  Where  we  have 
omitted  proofs,  the  reader  should  consult  Bishop’s  book,  since  his  proofs  are  usually  easily  adapted  to  our 
approach. 
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2.  The  Lambda  Calculus 


In  this  chapter  we  develop  the  syntax  and  semantics  of  the  lambda  calculus.  The  lambda  calculus  is  a 
model  of  computation  equivalent  in  power  to  the  Turing  Machine  and  other  well-known  models  of  compu¬ 
tation.  By  the  generally  accepted  definition  of  computability,  the  lambda  calculus  can  compute  any  com¬ 
putable  function. 

3.1  Syntax 

£.11  Token*: 

We  assume  a  denumerable  infinity  v  of  symbols: 

v  -  {a  ,b  ,e  ,  .  .  .  ,  z  ,aa,ab,  A  ,B ,  ,  succw  ,predw 

These  tokens  can  be  considered  strings  made  of  symbols  chosen  from  an  underlying  finite  alphabet  of  sym¬ 
bols. 

A  token  from  u  is  a  formula  of  the  lambda  calculus. 

Comment*:  Tokens  are  used  to  denote  constants  and  for  the  names  of  functions  and  operations  and  for 
the  formal  parameters  of  functions. 

£.1.2  Abstraction*: 

An  expression  of  the  form  (XvE),  in  which  v  is  a  token  and  £  is  a  formula  of  the  lambda  calculus,  is 
a  formula  of  the  lambda  calculus. 

Comment*:  The  intended  meaning  of  (A i/E)  is  the  function  that  takes  any  value  u  into  E .  For  exam¬ 
ple,  (As  (x  +  1))  is  the  successor  function,  since  it  takes  any  z  into  :  t  1. 

2.1.8  Application*: 

An  expression  of  the  form  (E  E  '),  in  which  E  and  £  '  are  formulas  of  the  lambda  calculus,  is  a  for¬ 
mula  of  the  lambda  calculus. 

Comment*:  The  notation  (/  z)  denotes  the  application  of  the  function  /  to  the  argument  z.  For 
example,  ((Az(z  +  1))  2)  denotes  the  application  of  the  successor  function  (A z  (z  -+■  1))  to  the  argument 
2. 
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8.1.4  Ezelusion: 


The  only  formulas  of  the  lambda  calculus  are  those  described  in  (1)  •  (3)  above. 

2.1.5  Definition  of  Bound  Identifier: 

A  variable  u  occurs  bound  in  a  formula  E  if  and  only  if  there  is  some  formula  F  such  that  (Xt/F)  is  a 
subformula  of  E . 

2.1.6  Definition  of  Free  Identifier: 

A  variable  occurs  free  in  a  formula  if  and  only  if  it  does  not  occur  bound  in  that  formula. 

3.3  Semantics 

2.2.1  Substitution  Notation: 

The  notation  ‘E  {v—A  }’  represents  the  result  of  substituting  the  formula  A  for  all  free  occurrences  of 
the  identifier  u  in  formula  E . 

2.2.2  Renaming  Rule: 

If  g  does  not  occur  free  in  E ,  then  we  can  perform  the  following  reduction  bg  renaming: 

(XzE)  =>  (Ay£{*-y}) 

2.2.3  Substitution  Rule: 

We  can  perform  the  following  reduction  bg  substitution: 

({XzE)  A)  -a  E{z~A  } 

provided  that  it  does  not  result  in  any  free  variable  of  A  becoming  bound. 

2.2.4  Definition  of  Reduction: 

A  reduction  is  a  sequence  of  one  or  more  applications  of  the  renaming  and  substitution  rules.  We 
write  E  =*  F  to  denote  that  there  is  a  reduction  from  E  to  F .  We  also  use  this  symbol  to  introduce 
abbreviations  and  explicit  definitions.  Note  that  this  means  that  all  abbreviations  and  definitions  are  in 
principle  eliminatable. 
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£.2.5  Definition  of  Normal  Reduction  Order: 

We  say  that  a  reduction  is  performed  in  normal  reduction  order  if  and  only  if  the  outermost  reducible 
application  is  always  reduced  before  any  inner  applications.  That  is,  we  never  reduce  a  formula  by  substi¬ 
tution  that  is  a  subformula  of  an  application  that  could  be  reduced  by  substitution. 

2.2.6  Definition  of  Normal  Form: 

A  formula  is  said  to  be  in  normal  form  if  and  only  if  the  substitution  rule  cannot  be  applied  to  it,  even 
after  intervening  applications  of  the  renaming  rule. 

2.2.7  Church- Rotter  Theorem: 

If  a  formula  has  a  normal  form,  then  that  normal  form  is  unique  (up  to  applications  of  the  renaming 
rule).  Thus,  the  normal  form  of  a  formula  can  be  considered  that  formula’s  value. 

2.2.8  Definition  of  a  Formula  Being  Defined: 

We  say  that  a  formula  is  defined  if  and  only  if  it  has  a  normal  form.  Note  that  by  Turing’s  Halting 
Theorem  it  is  not  finitely  decidable  (i.e.,  algorithmically  decidable)  whether  or  not  a  formula  is  defined. 

2.2.9  Univertality  of  Normal  Reduction  Order: 

If  a  formula  is  defined,  then  a  normal  order  reduction  will  reach  its  normal  form. 

2.2.10  Definition  of  Convertibility: 

We  call  two  formulas  E  and  F  convertible ,  and  write  'E  w*  F\  if  and  only  if  there  is  a  normal  form 
N  such  that  E  =*  N  and  F  =*  N .  That  is,  two  formulas  are  convertible  if  and  only  if  they  have  the 
same  normal  forms. 

Note  that  if  E  and  F  are  defined  then  E  F  is  finitely  decidable:  simply  reduce  both  E  and  F  to 
their  normal  forms  (which  are  finite  formulas),  and  compare  them.  This  comparison  can  be  done  finitely 
since  the  formulas  are  finite. 
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3.  Basic  Definitions 


In  this  chapter  we  introduce  a  number  of  abbreviations  and  notations,  which  make  the  lambda  cal¬ 
culus  convenient  to  use  as  a  functional  programming  language.  These  extensions  include  one  atomic  data 
type,  Booleans,  and  one  composite  data  type,  LISP-style  lists. 

5.1  Parameters 

3.1.1  Dyadic  Function:  (A(xy)£)  =»  (Xz(XyE)) 

3.1.2  \-adic  Function:  (A(x,x2  ••■*»)£)  =r>  (Azt(Azj(  -  •  (A x„£)  •••))) 

3.1.3  Omittion  of  Parentheses  (Application): 

Dyadic:  (/  x  y)  =»  ((/  x)  y) 

N-adic:  (/  x,  x,  z. )  =*  (  ((/  *i)  *i)  ‘  *» ) 

3.1.4  Omittion  of  Parenthetet  (Abitr action): 

XzE  =*  (Ax£ ) 

A(x„  .  .  .  ,  x„)£  =*  (A(x,,  .  .  .  ,  x„)£) 

3.1.5  Alternative  Brackett:  Other  forms  of  parentheses,  such  as  [~]  and  {“},  may  be  used  in  place  of  (— ) 
for  readability. 

5.1  Com  bins  tors 

3.2.1  Identity:  Id  »o  Axx 

CommenU:  Note  that  Id  a  a  ,  for  all  a  . 

3.2.2  Contant:  K  •*  A (kz)k 

Comment*:  The  function  ‘K  e  ’  returns  the  value  e  for  all  arguments: 

K  e  a  =*  X(kz)k  c  a  =*  c 

3.2.3  Projection  Function: 

1st  ^  A(xy)x 
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2nd  =•  A(ry)y 


Comment*:  Thus,  1st  a  b  ^  a  and  2nd  a  b  =»  b .  These  functions  are  usually  used  in  conjunc¬ 
tion  with  other  combinators. 

3.2  4  Monadic  Competitor: 

B  =*  A (fgz)\f(gz)\ 
f  '9  =*  B  f  9 

Comment*:  Both  'B’  (the  prefix  form)  and  ‘  *  ’  (the  infix  form)  compute  the  composition  of  two 
monadic  functions.  To  see  this  note  that 

(log*sin)a  ^  B  log  sin  a  log  (sin  a  ) 

3.2.5  Monadic  Formalizer: 

*  =*  A(/a6*)j/  (oi)(6r)] 

/  *  (“.*»)  =*  **[/  («*)(**)! 

Comment*:  If  ‘sum’  is  the  addition  function,  then  *sum‘(/  ,g)'  represents  the  functional  sum  of  / 
and  g ,  since 

sum  ’  (f  ,g)  *  ^  sum  (f  *)  {g  *) 

For  example,  sum  "  (sin, cos)  is  a  function  to  compute  the  sum  of  the  sines  and  cosines  of  its  argument: 

sum  *  (sin, cos)  9  sum  (sin  9  )  (cos  9  ) 

Thus  this  operation  feeds  a  single  input  through  two  monadic  functions,  and  then  feeds  the  outputs  of 
these  into  a  dyadic  function. 

3.2.6  Dyadic  Competitor: 

T  =*  A(/g){A(*y  )(/  (pry)!) 
f  9  -  r  f  9 

Comment*:  For  example, 

(log:sum)  a  I  •*  X(zy  )(log(sum  r  y)j  a  b  log(sum  a  b) 
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This  composes  a  monadic  function  with  the  output  of  a  dyadic  function. 

3  2  7  Dyadic  Formalizer: 

T  =»  X(fabzy)\f  (azy)(bzy  )| 

f  (a,b)  =*  A(zy )/  (axy)(6ry)| 

Comment*:  This  operation  feeds  two  arguments  to  each  of  two  dyadic  functions,  and  then  feeds  the 
two  results  to  a  third  dyadic  function.  Thus,  if  /  and  g  are  dyadic  functions,  then  sum:(/  g )  is  the 
functional  sum  of  /  and  g.  For  example,  to  compute  (z-y)(z-i-y)  we  can  use  the  function 
prod:  (dif, sum): 

prod:(dif,sum)  z  y  «■*  prod  (dif  z  y)  (sum  z  y) 

3.2.8  Dual  Competitor:  fi  =*  \{fghzy  )[f  (gz  )(hy  )] 

Comment*:  This  preprocesses  the  two  inputs  to  a  dyadic  function  with  two  separate  monadic  func¬ 
tions.  For  example, 

(fi  sum  sin  cos)  a  b  =»  A(/yAzy)[/  («*)(*?)!  sum  sin  cos  a  6  =■»  sum  (sin  a  )  (cos  b  ) 

Thus,  ‘fl  sum  sin  cos’  can  be  read  “the  sum  of  the  sin  of  the  first  and  the  cosine  of  the  second.’’  This  com¬ 
bination  can  also  be  expressed  by  the  dyadic  compositor  and  the  projection  functions: 

sum:(sin  *  1st,  cos  *  2nd) 

To  see  this,  observe 

sum: (sin  *  1st,  cos  •  2nd)  a  6  sum  (sin  *  1st  o  )  (cos  •  2nd  a  b  ) 

=*  sum  (sin[lst  a  6])  (cos  2nd  a  b  j) 

-*  sum  (sin  a  )  (sin  b  ) 

3.2.9  Elementary  Permutator  (Convection):  conv  ^  A ( /  xy  ) ( / yx  ) 

Comment*:  For  example,  if  (dif  z  y)  computes  the  difference  of  z  and  y ,  then  (conv  dif  1)  is  the 
predecessor  function,  since 


(conv  dif  1)  a  ”*  A(/zy  )(/yz)  dif  1  a  ^  dif  a  1 


Thus  conv  exchanges  the  arguments  of  a  dyadic  function. 

3. £.10  S  Combinator:  S  =*  X(fgz)\fz(gz)  j 

Comment*:  The  S  combinator  is  difficult  to  explain.  It  is  best  understood  as  taking  a  single  argument, 
preprocessing  it  by  a  monadic  function,  and  then  passing  both  the  preprocessed  and  unpreprocessed  ver¬ 
sions  to  a  dyadic  function.  For  example,  (S  prod  log)  multiplies  a  number  by  its  logarithm: 

(S  prod  log)  p  =*  prod  p  (log  p  ) 

3.2.11  Pit  Combmator:  *  =»  X(fgzy  )[/  (gz  )(gy  )) 

Comment*:  This  combinator  preprocesses  both  arguments  of  a  dyadic  function  by  the  same  monadic 
function.  Thus  (It  sum  log)  is  a  function  to  add  logarithmically: 

(4  sum  log)  a  b  =>  sum  (log  a  )  (log  b  ) 

Note  that 

♦  /  g  <*=>  A  f  g  g  <=*>  /  :(g  ‘1st,  g  ‘2nd). 

3.2.12  Paradoxical  Combinator:  Y  =*  Xf  {  Xz\f  (xx)J  A  *  [  /  (zz )]  } 

Comment*:  The  Y  combinator  is  difficult  to  explain.  This  definition  of  Y  satisfies  the  functional  equa¬ 
tion: 

YF  ^  F(Y  F)  F(F(YF))  -a  •  •• 

To  see  this  observe 

YF  Ax(F(xx)|  Ax(F(xx)|  -•  F{Ax[F(xx)  Ax[F(xx)j}  ••• 

Thus  it  can  be  considered  the  fized-point  finding  operation. 

3.2.13  Recurtive  Definition:  *rec  F  =  D  ’  ‘F  ^  Y  A FD  ’ 

Comment*:  For  example,  a  definition  of  the  form 

,  /n  =  0  -  1 

rec  fac  =  An)  ,  ,,  , 

\*1**  «  x  |fac(n  -  l)j 

is  to  be  replaced  by 
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f&c 


(  n  -0  *1 

Y  Afac  An^eUe  n  x  f&c  _  ,j 


The  Y  combinator  in  effect  embeds  the  definition  of  fac  in  itself. 

3.8.14  Seminormal  Form:  We  say  that  a  formula  is  in  teminormal  form  if  and  only  if  the  only  reductions 
that  can  be  performed  on  it  and  its  descendants  are  of  the  form  Y  F  =*  F  (  Y  F). 

Comment*:  Because  of  their  use  of  the  paradoxical  combinator  many  useful  formulas,  for  example  the 
recursive  definition  of  ‘fac’  exhibited  above,  do  not  have  a  normal  form.  Conceptually,  however,  this  for¬ 
mula  represents  a  value,  viz.,  the  factorial  function.  If  we  view  any  application  of  Y  as  being  a  potential 
application,  that  is  only  performed  when  necessary  to  allow  the  reduction  of  non-Y  applications,  then  we 
can  see  that  a  formula  is  in  seminormal  form  when  none  of  its  actual  applications  can  be  reduced.  That 
is,  a  formula  is  in  seminormal  form  if  either: 

1.  it  is  in  normal  form,  or 

2.  it  is  not  in  normal  form,  but  the  only  applications  that  can  be  reduced  are  Y  applications,  and  per¬ 
forming  these  Y  applications  would  not  enable  any  non- Y  applications. 

Thus  a  formula  in  seminormal  form  can  be  thought  of  as  dormant.  It  can  be  activated  by  using  it  in  an 
application.  For  example,  if  A  (zy)£  is  in  normal  form,  then  Y[A(zy  )E  ]  is  in  seminormal  form.  This 
function  can  be  activated  by  applying  it  to  an  argument  a ,  since  this  application  allows  non-Y  reductions 
to  take  place: 

{Y(A(zy )E |}  a  =*  A(zy )E  {Y[A(zy )£J}  a 

=•>  E{z  -  {Y(A(zy)£j},  y  -  □} 

We  loosely  use  the  term  defined  to  mean  that  a  formula  has  either  a  normal  or  seminormal  form,  and  only 
discriminate  between  the  two  when  necessary. 

3.3  Conditionals 

8.3.1  True  Value:  true  =*  A(zy)z 
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Comment*:  Note  that 


true  a  b  =*  A(zy)z  a  b  ^  a 
Thus  true  selects  the  Tint  of  Us  two  arguments. 

3.3.2  False  Value:  false  A(zy)y 
Comments:  Note  that 

false  a  6  -*  A(zy)y  a  b  =*  6 
Thus  false  selects  the  second  of  its  two  arguments. 

3.3.3  Definition  of  Boolean  Value: 

1.  The  normal  form  of  ‘true'  is  a  Boolean  value. 

2.  The  normal  form  of  ‘false’  is  a  Boolean  value. 

3.  The  only  Boolean  values  are  those  described  in  (1)  and  (2)  above. 

3.3.4  Negation:  not  Az(z  false  true) 

Comments:  For  example, 

not  true  =■*  Az(z  false  true)  true  true  false  true  ^  A(zy  )z  false  true  -*  false 

3.3.5  Conjunction:  and  =*  A(zy)(z  y  false) 

Comments;  For  example, 

and  true  false  «•  true  false  false  •*  A(zy)z  false  false  -•  false 
3.3.0  Disjunction:  or  =■*  A(zy  )(z  true  y ) 

Comments:  For  example, 

or  false  true  =*>  false  false  true  ^  true 
3.3.7  Equivalence:  equiv  =*  S  z  not  y 
Comments:  For  example, 

equiv  true  true  S  true  not  true  «*  true  true  (not  true)  =*  true 

equiv  false  true  -*  S  false  not  true  -*  false  true  (not  true)  ^  not  true  -*  false 
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3.3.8  Dyadic  Conditional: 


P  ~C 

else  A 


P  C  A 


Comments:  This  depends  on  the  fact  that  if  P  is  a  Boolean  value,  then  it  will  choose  between  C  and 
A  depending  on  whether  it  is  true  or  false. 


3.3.9  ,Y-  adie  Conditional: 

Pi-Ct 

P^C, 

=*  p,  <M j 

Pn^CK 

else  A 

3.4  Lists 

3.4.1  First  of  Three:  lofS  =>  A (zyz)z 
Comments:  Note  that  for  any  a  ,  b  ,  e  : 


Pi  Ci  (P,  c,  •  (P»  Cn  A)  ■  •  ) 


lof3  a  b  e  =>  [A(ry*)ri  a  6  e  =»  a 

Notice  that  lofS  selects  the  first  of  three  values  in  much  the  same  way  that  true  selects  the  first  of  two 
values.  In  effect  the  formulas  lof3,  2of3  and  SofS  form  the  truth  values  of  a  three-valued  logic. 


3.4-S  Second  of  Three:  2of3  =»  A (zyz)y 
3-4-3  Third  of  Three:  3of3  =»  A(zyz  )z 
3-4-4  Value:  list  =*  A(nht)[A«(s  n  h  t )] 

Comments:  This  should  be  thought  of  as  a  triple  (n  ,h  ,t ),  in  which  n  is  the  null-flag  and  h  and  t  are 
the  head  and  tail  of  the  list.  The  latter  have  significant  values  only  if  the  null-fiag  is  false.  The  parame¬ 
ter  s  can  be  thought  of  as  a  potential  selector,  which  will  select  one  of  n  ,  h  or  t .  Notice  that 

list  false  A  B  =*  A#  («  false  A  B) 


•  »  •  .  •  „  * .  1  % 

V  %*  '  *  '  .•  A  n*  -*  V  . 


.v  .*• 


3-4-5  \ull  Lilt:  nil  =*  list  true  j_  _ 

Comment i:  The  first  component  of  the  triple,  true,  indicates  that  the  list  is  null,  and  that  the  remain¬ 
ing  two  components  are  not  used.  The  symbol  ‘j_’  has  no  special  significance,  although  it  is  intended  to 
suggest  an  undefined  value;  it  is  merely  a  token  like  any  other  in  the  lambda  calculus.  Note  that 

nil  =»  list  true  j_  j_  \i  (i  true  j_  j_) 

which  is  the  normal  form  for  the  null  list. 

3  4-6  Lilt  Construction;  cons  A (ht ) (list  false  h  t ) 

Commenti:  For  example, 

cons  A  B  =*  list  false  .4  B  =»  As  (s  false  A  B) 

That  is,  the  result  of  ‘cons  A  B’  is  a  list  triple,  whose  first  element  is  false  (meaning  the  list  is  nonnull), 
and  whose  second  and  third  elements  are  A  and  B ,  the  head  and  tail  of  the  list. 

3.4-7  Suit  Teit:  null  =>  Az(z  lof3) 

Commenti:  This  simply  selects  the  first  (null-flag)  component  of  a  list  value; 

null  nil  =■*  nil  lofS  =*  As  (s  true  j _ l)  lofS  =*  lofS  true  j_  jl  <=»  true 

3-4  3  Head  of  Lilt:  hd  =*  Ax  (z  2ofS) 

Commenti.  This  simply  selects  the  second  of  the  three  components  of  a  list  value; 
hd(cons  A  B )  =>  (cons  A  B )  2of3  =*  As  (s  false  A  B )  2ofS  2of8  false  A  B  =*  A 
Note  that  hd  nil  «=*  j_. 

3-4-9  Tail  of  Lilt:  tl  =»  Az(z  3of3) 

Comment!.  This  simply  selects  the  third  of  the  three  elements  of  a  list  value.  We  have  that 
tl(cons  A  B  )  B  .  Similarly,  it  is  easy  to  show  that  if  L  is  a  nonnull  list  then 

cons  (hd  L )  (tl  L  )  «=*  L 

Simply  observe  that  if  L  is  a  nonnull  list  then  it  has  the  form  (list  false  A  B ).  Note  also  that 
tl  nil  j_. 
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3.4- 10  Lilt  Notation: 

<ii,  *j,  •  •  >  =»  cons  i!  (cons  z2  ■  •  •  (cons  z„  nil)  •  •  ) 

More  formally, 

<  >  =*>  nil 

<zt,  r j,  ..,*«  >  =»  cons  zt  <z2,  .  .  .  ,  z„  > 

Notice  that  by  previously  derived  properties: 

hd  <*,,  z2,  .  .  .  ,  *B  >  =*  hd(cons  z,  <z2,  .  .  .  ,  zn  >)  =»  z, 

tl  <Zj,  Zj,  .  .  .  ,  z„  >  =»  tl(cons  Zj  <z2,  .  .  .  ,  z„  >)  =»  <z2,  .  .  .  ,  z„  > 

3.4- 11  Pair  Lift:  pair  =>  A(zy)<z,y> 

Comment*:  The  formula  ‘pair  z  y  ’  constructs  a  two-element  list  <z  ,y  >.  Note  that  ‘pair  a ’  is  a  func¬ 
tion  that  constructs  a  list  whose  first  element  is  a  : 

(pair  a)  y  =*  pair  ay  =>  <0  ,y  > 

3.4  12  Definition  of  List: 

1.  The  normal  form  of  ‘nil’  is  a  lilt. 

2.  If  L  is  a  list  and  z  is  any  formula  having  a  normal  form,  then  the  normal  form  of  ‘cons  z  L  ’  is  a 
lilt. 


3.  The  only  lists  are  those  described  in  (1)  and  (2)  above. 

S.4  IS  Example  Normal  Form: 

For  concreteness  we  exhibit  the  normal  form  of  a  particular  list.  Suppose  ‘1’  and  ‘a '  are  two  tokens 
(symbols  in  v).  Then  we  derive  the  normal  form  of  the  list  ‘<l,o  >’  as  follows: 

<l,o  >  =*  cons  1  (cons  a  nil) 

=»*  (list  false  1  (cons  a  nil)) 

(list  false  1  (list  false  o  nil)) 

(list  false  1  (list  false  a  (list  true  j_  _))) 
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1 


I 


»*  A#  (»  false  1  (list  false  a  (list  true  j _ l))) 

=*  A*  (s  false  1  A*  (s  false  a  (list  true  j_  _l))) 

=»  A*  («  false  1  As  (s  false  a  As  (s  true  j _ >_))) 

=»  As  (s  A(xy)y  1  As  (s  A(xy)y  a  As  (s  A(*y  )x  j _ ))) 

The  final  formula  above  is  strictly  speaking  not  in  normal  form,  which  would  require  replacing  the 
suppressed  parentheses.  We  omit  this  step  as  it  would  make  the  formula  too  unreadable. 

3.4.14  Nonstandard  Objects: 

This  representation  permits  several  kinds  of  nonstandard  lists.  In  particular,  the  structure 
(list  true  x  y }  behaves  exactly  like  nil,  the  empty  list,  but  contains  two  pieces  of  information,  namely  x 
and  y  .  Also,  although  in  the  application  (cons  x  y )  the  argument  y  is  expected  to  be  a  list,  this  is  not 
necessary.  Whatever  it  is,  the  tl  function  will  extract  it. 


I 

i 


I 
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4.  Natural  Numbera 


In  this  chapter  we  define  a  representation  for  natural  numbers  and  for  the  basic  operations  on  natural 
numbers.  Natural  numbers  are  represented  in  unary  notation  as  lists  of  identical  tokens. 

4.1  Primitive  Ideas 

4-1.1  Definition  of  Natural  Number: 

A  natural  number  is  a  (possibly  null)  list,  all  of  whose  elements  are  nil. 

4-1.2  Numeric  Denotation*: 

0  =>  <  > 

1  =»  <nil> 

2  =*>  <nil,nil> 

3  =s»  <nil,nil,nil> 

4  =>  <nil,nil,nil,nil> 

9  =*  <nil,nil,nil,nil,nil,nil,nil,nil,nil> 

4-1.3  Example  Normal  Form: 

For  concreteness  we  exhibit  the  normal  form  of  a  particular  natural  number: 

0  =»  nil 

=®  (list  true  j_  j_) 

=>  A«  («  true  j_  j_) 

=*  A*  (#  A(zy  )x  j_  j_) 

4-1-4  Nonstandard  Object*: 

There  are  no  nonstandard  natural  numbers,  since  they  are  defined  to  be  lists  of  nils.  However,  it  will 
be  seen  that  the  arithmetic  operations  and  relations  depend  only  on  the  lengths  of  these  lists,  so  that  lists 
of  things  other  than  nils  can  be  considered  nonstandard  natural  numbers.  For  example,  the  normal  forms 
of  <_>,  <tme>,  < false > ,  <<nil>>,  etc.  could  be  considered  alternate,  nonstandard  representations  of 


I. 


4-1.5  Zero  Predicate:  zero? y  =>  null 

Comment s:  For  example, 

zero ?  v  0  =s>  null  <>  =»  true 
zero?y  1  =*  null  <nil>  =>  false 

4-1.6  Successor:  succ y  =s>  cons  nil 

Comments:  Since  the  natural  number  n  is  represented  by  a  list  of  n  nils,  the  successor  function  simply 
adds  another  nil  to  the  list.  For  example, 

succy  3  =>  cons  nil  <nil,nil,nil>  <nil,nil,nil,nil>  «=>  4 

^.1.7  Predecessor:  predy  ==»  tl 

Comments:  For  example, 

predy  4  =*  tl  <nil,nil,nil,nil>  =>  <nil,nil,nil>  <=>  3 

4.3  Arithmetic 

4-2.1  Repeated  Composition: 

{zero?y  n  -*  a 
else  /  [rpt  /  (predy  n  )  a  ] 

Comments:  Notice  that: 

rpt  /  0  z  =>  z 

rpt  /  1  z  =*  /z 

rpt  /  2  z  =*  /  ( /  x  ) 

rpt  /  3  z  =»  /  (/  (/z  )) 

rpt  /  n  z  /  "  z 

In  general,  rpt  /  n  is  the  n  -fold  composition  of  /  with  itself. 
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4-2.2  Sum:  sum#  =*  rpt  succ # 

Comments:  Thus,  the  sum  of  m  and  n  is  the  m  th  successor  of  n  : 

sum#  m  n  =9  rpt  succ#  m  n  =»  succ#mn  =>  m-i-n 
4-2.3  Difference:  dif#  =*  conv(rpt  pred# ) 

Comments:  The  effect  of  this  definition  is  that  the  difference  of  m  and  r»  is  the  n  th  predecessor  of  m  : 

dif#  m  n  =9  conv (rpt  pred# )  m  n  =*  rpt  pred#  n  m  =»  pred#*m  =»  m  -  n 

4-2-4  Product:  prod#  =*  A ( mn  ) [rpt  (sum#  m  )  n  Oj 

Comments:  Notice  that  (sum#  m  )  is  the  function  that  adds  m  .  Hence  the  product  of  m  and  n  is 
defined  to  be  the  result  of  adding  m  to  zero  n  times: 

prod#  m  n  =>  rpt  (sum#  m )  n  0  =*  (sum#  m)*0  =*  m-t-m+  ■  •  •  +m+0  =»  mn 

4.3  Predicates  and  Relations 

4-3.1  Less- Than  Relation: 

zero?#  n  -*  false 
rec  It#  =  A(m«  )  zero?#  m  -*  true 

else  It#  (pred#  m)  (pred#  n) 

Comments:  This  definition  is  based  on  the  fact  that 

m  <  n  <*>  m  - 1  <  n  - 1 

Thus,  It#  m  n  operates  recursively,  subtracting  one  from  each  of  m  and  n  until  one  or  the  other  reaches 
0,  at  which  point  the  truth  value  of  the  relation  is  known: 

It#  m  0  =*•  false 

It#  0  n  =*  true,  if  n  nonzero 

It#  m  n  =»  It#  (m-l)(n-l),  if  m ,  n  nonzero 

4-3.2  Greater- Than  Relation:  gt#  =*  conv  It# 

Commente:  Observe  that 
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gt|v  m  n  =*  conv  Itjy  n  r  =»  Itjv  n  m 
4- 3. 3  Inequality  Relation:  tie#  =»  or:(ltjv  ,  gt/v  ) 

Comments:  Two  natural  numbers  are  unequal  if  one  is  either  greater  or  less  than  the  other: 

nejv  m  n  =*>  or:(lt/y  ,  gt# )  m  n  *9  or(lt^  m  n )  (gtw  m  n) 

4-3-4  Equality  Relation:  eqN  =*  not:ne  v 

Comments:  Two  natural  numbers  are  equal  if  they  are  not  unequal: 

eqjy  m  n  =9  not:nelV  m  n  =9  not(new  m  n) 

4-3.5  Less-or-Equal  Relation:  lejy  =»  not:gtN 

Comments:  One  number  is  less  than  or  equal  to  another  if  it  is  not  greater  than  the  other: 

lejy  m  n  =»  not:gt^  m  n  =>  not(gt,v  m  n ) 

4-3.6  Greater- or-  Equal  Relation:  geN  =*  not: It N 
4.4  Miscellaneous  Operations 

{ltw  m  n  —  0 
elsesucc^quo*  (dif*  m  n)  n 

Comments:  The  definition  of  quotient  is  baaed  on  the  following  identity: 

m-T-n  =  1  +  (m  —  n  )  -5-  n 


Note  that  if  m  is  greater  or  equal  to  n  then 


quow  m  N  =9  succjyiquojy  (dif*  m  n )  n  =*>  succyjquojv  dif^v  m  n)  nj 


(zero?*  n  -»  1 
else  prod*  {n  ,  [fac  (pred*  n)\} 


4  -4-3  Definition  of  Sequence: 

A  sequence  is  any  formula  N  such  that  for  each  nonzero  natural  number  k  the  formula  ‘Nk  ’  is 
defined. 
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If  C  is  any  class  of  normal  forms  and  N  is  a  sequence,  then  N  is  called  a  tequenee  of  C*  if  the  normal 
form  of  each  formula  W  is  in  C  . 

4-4-4  Subieript  Notation: 

If  .V  is  a  sequence  we  usually  write  Nk  for  Nk : 


S.  Integer* 


In  this  chapter  we  define  integers  as  pairs  comprising  a  sign  (represented  by  a  Boolean  value)  and  a 
magnitude  (represented  by  a  natural  number). 

5.1  Primitive  Ideas 

5.1.1  Definition  of  Integer: 

An  integer  is  a  pair  (two  element  list),  whose  first  element  is  a  Boolean  value  and  whose  second  ele¬ 
ment  is  a  natural  number. 

5.1.2  Attach  Plut-tign:  plus/  =>  pair  false 

Comment*:  This  operation  converts  a  natural  number  into  a  plus-signed  integer  by  attaching  a  plus- 
flag  (false  value)  to  it.  For  example, 

plus/  S  =*  pair  false  3  =*■  < false, 3 > 

5.1.3  Attach  Minu*-*ign:  minus/  =*  pair  true 

Comment*:  This  operation  converts  a  natural  number  into  a  minus-signed  integer. 

5.1-4  Denotation*: 

-*-0  =>  plus/  0 
—  1  =*  plus/  1 
-2  =»  plus/  2 
-1  =>  minus/  1 

5.1.5  Example  Normal  Form: 

For  concreteness  we  exhibit  the  normal  form  of  ‘+0’: 

*0  =*  plus/  0 
=»  <  false,  0> 

=*  <A(ty)y,As(#  A(zy)z  il)> 

=*  As  (s  A(zy)y  A(zy)y  As  (s  A(zy)y  As  (s  A(zy)z  j_  Jj  As  (s  A(zy  )z  _l))) 
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S.J.6  Nonstandard  Objects: 


This  representation  permits  one  nonstandard  object,  -0.  This  is  because  ->-0  and  -0  have  different 
representations,  <false,0>  and  <true,0>,  although  the  arithmetic  operations  and  relations  treat  them  as 
equivalent. 

5.1.7  Extract  Magnitude:  mag;  =*  hd  •  tl 

Comments:  This  operation  converts  an  integer  into  a  natural  number  by  discarding  its  sign:  For  exam¬ 
ple, 

mag;  +1  =*  mag;  < false,  1>  =*  hd  •  tl  <false,l>  =*  hd(tl  <false,l>)  =>  hd  <1>  =*  1 
mag;  -1  =»  mag;  <true,l>  =*  hd  •  tl  <trae,I>  =»  hd(tl  <trne,l>)  =>  hd  <1>  =*  1 

S.2  Primitive  Predicates 

5.2.1  Minus-sign  Test:  minus?;  =*  hd 

Comments:  This  operation  tests  if  an  integer  bears  a  minus  sign: 

minus?2  +1  =*  minus?;  <false,l>  =*  hd  <false,l>  =>  false 

Notice  that  a  number  bearing  a  minus  sign  is  not  the  same  thing  as  it  being  negative,  since  -0  bears  a 
minus  sign.  True  positive  and  negative  tests  are  defined  later. 

5.2.2  Plus-sign  Test:  plus?;  =>  not  "minus?; 

Comments:  For  example, 

plus?;  +1  =*  not  "minus?;  +1  =*  not(minus?;  +1)  =»  not  false  =»  true 
Again,  bearing  a  plus  sign  is  not  the  same  as  being  positive,  since  +0  bears  a  plus  sign. 

5.2.3  Zero  Test:  zero?;  =*  zero?*  ’mag; 

Comments:  For  example, 

zero?;  +1  =*  zero 7 N  "mag;  +1  zero?*  (mag;  +1)  =»  zero?*  1  =*  false 
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1.3  Primitive  Operations 


5.3  1  Absolutt  Value:  aba*  ^  plus/  •  mag/ 

Comment*:  For  example, 

abs/  -1  =*  plus/ (mag/  -1)  =»  plus/  1  <*=*•  ->-l 

5.S.S  Segate:  neg/  pair*  (plus?/  , mag/ ) 

Comment*:  This  operation  negates  an  integer.  For  example, 

neg/  -1  =*  pair*  (plus?/ ,mag/)  -1 

=>  pair  (plus?/  -1)  (mag;  -1)  pair  false  1  <=>  +1 


5.3.3  Sueee**or:  succ/ 


zero: 2  n 
Arc  plus? 2  n 
minus?/  n 


-  +1 

-»  plus/  *su ccw  'mag2  n 
-*  minus/  *  pred  y  *  mag/  n 


Comment*:  This  operation  computes  the  successor  of  an  integer.  The  main  complication  is  to  avoid 
treating  -0  and  +0  as  different  integers.  The  cases  are: 

succ/  ±0  =*  +1 

succ/  +n  =>  plus/  [succ y  (mag/  +  n  )|  =*  +  (n+l) 
succ/  -n  =>  minus/ [predy  (mag/  — n)|  =»  — (n-1) 

5.3  4  Predeeeuor:  pred/  =»  neg/ *succ/ *  neg/ 

Comment*:  Notice  that 

pred/  n  =»*  neg/ *  succ /*  neg/  n  =*  neg/ (succ/ (neg/  n)j 
=  -[1  +  (-  n  )]  =  — [1— r»  ]  =  n  -1 
/ 

sero?/  n  -*  m 

5.3.5  Sum:  rec  sum/  =  A(mn)'  minus?/  n  -»  pred/  [sum/  m  (succ /  n)j 

plus?/  n  -♦  succ/  [sum/  m  (pred/  r» )) 


Comment*:  The  sum  of  integers  is  defined  recursively  by  the  following  cases. 
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m  0  =  m 

m+n  =  [m  +  («  +  l)j  -  1,  if  n  <0 
m  +  n  =  [m  +  (»  -1)|  1,  if  n  >0 

5.3.6  Difference:  dif/  =»  sum*  : ( 1st,  neg/  •  2nd) 

Comment i:  Note  that 

dif/  m  n  =*  sum/ :( 1st,  neg/  •  2nd)  m  n  =*  sum/  m  (neg/  n) 

=>  sum/  m  (neg/  n)  =  m+(-n)  =  m-n 

S.4  Predicates  and  Relations 

5.4- 1  Equality:  eq/  =»  zero? /  :dif/ 

Comment*.  This  depends  on  the  fact  that  m  =  n  if  and  only  if  m  -  n  is  zero: 

eq/  m  n  =*  zero?/:dif/  m  n  ^  zero?/(dif/  m  n) 

5.4.2  Negative  Predicate:  neg?/  =0  and  “  (minus?/  ,  not  •  zero?/ ) 

Comments:  This  is  a  true  negative  test  (in  contrast  to  minus?/ ),  and  depends  on  the  fact  that  m  <0  if 
and  only  if  m  bears  a  minus  sign  and  is  nonzero: 

neg?/  m  =>  and  •  (minus?/  ,  not  •  zero?/ )  m  =*  and  (minus?/  m  )  jnot(zero?/  m  )j 

5.4- 3  Positive  Predicate:  pos?/  =0  andserie  (plus?/ ,  not  *  zero?/ ) 

Comments:  An  integer  is  positive  if  it  is  nonzero  and  bears  a  plus  sign: 

pos?/  m  =*  and  •  (plus? / ,  not  •  zero? / )  m  =»  and  (plus?/  m )  [not(zero?/  m  )] 

5.4- 4  Less- Than  Relation:  It/  =>  neg?/:dif/ 

Comments:  We  make  use  of  the  fact  that  m  <n  if  and  only  if  m  -n  is  negative: 

It/  m  n  =*  neg?/:dif/  m  n  =*  neg?/(dif/  m  n ) 

5.4- 5  Greater- Than  Relation:  gt /  =*  pos?/:dif/ 

Comments:  Note 
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gt/  m  n  =»  pos?/.dif/  m  n  **  pos? /  (dif/  m  n) 


5-4-d  Greater- or-  Equal  Relation:  ge /  =»  not: It/ 


Cofnm<nt«.  Here  we  make  use  of  the  fact  that  m^n  is  true  if  and  only  if  m  <n  is  false: 


ge/  m  n  =»  not:lt/  m  n  =*  not(lt2  m  n) 


5.^.7  Leu- or- Equal  Relation:  \ez  =*  not:gt/ 


5.4  S  Maximum:  max/  =»  A(ry) 


/  (g *Z  * 

^else  y 


y)  -  * 


S.S  Multiplicative  Operation! 


5.5.1  Product: 


j equiv  (neg/  m)  (neg/  n)  —  plus/  .prod N  (mag/  m)  (mag/  n) 
prod/  =*  '^(,nn  J^else  minus/ :prody  (mag/  m  )  (mag/  n) 

Comments:  This  definition  operates  by  attaching  the  appropriate  sign  to  the  result  of  a  natural 
number  multiplication  of  the  magnitudes  of  the  factors.  The  only  complication  is  to  take  care  of  the  sign 
of  the  result.  In  this  connection  note  that 

equiv  (neg/  m )  (neg/  n ) 

is  true  if  and  only  if  the  signs  of  m  and  n  are  the  same. 


5.5.2  Quotient: 


A(mn ) 


(equiv  (1 

else  mi 


'  (neg/  m)  (neg/  n)  -*  plus/:quow  (mag/  m)  (mag/  n) 
minus/  :qu o/f  (mag/  m  )  (mag/  n  ) 


Comment*:  Like  the  product  operation,  quotient  operates  by  attaching  the  appropriate  sign  to  the 


result  of  the  corresponding  natural  number  operation. 
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6.  Rational  Number* 


In  this  chapter  rational  numbers  and  their  operations  are  defined,  based  on  the  representation  of  a 
rational  number  as  a  pair  of  integers. 

6.1  Primitive  Idea* 

6.1.1  Definition  of  Rational  Number: 

A  rational  number  is  a  pair  (two  element  list),  each  of  whose  elements  is  an  integer. 

6.1.2  Conttruction:  rat<j  =>•  pair 
Comment*:  Note  that 

ratg  m  n  =?>  pair  m  n  =>  <m,n> 

6.1.3  Notation:  m  /  n  =>  <m  ,n  > 

6.1.4  Denotation*: 

m  /  1  =>  in  /  +1 
0/  1  =»  +0/  1 

1/  1  =>  +1/1 

2/  1  =>  +2/  1 

1/  2  =•>  +1/  +2 

6.1.5  Example  Normal  Form: 

For  concreteness  we  exhibit  the  normal  form  of  0/0: 

0/  0  =*  <+0,+0> 

=*  As  (*  false  +0  A*  (1  false  *0  As  (1  true  j_  j_ ))) 

=*  A#  (*  A(zy)y  A*  (t  A(zy  )y  A(zy)y  A*  (s  A(zy)y  A*  (#  A(zy)z  j _ i_) 

A#  (s  A(zy)z  J_  J_))) 

A*  (s  A(zy)y  A#  (*  A (zy)y  A(zy)y  A*  (*  A(zy)y  A t  (*  A(zy)z  j_  J_) 

A#  (s  A(zy  )z  jl  J_)))  A#  (s  A(zy  )z  J_  _]_))) 

The  normal  form  of  1/2  is  shown  in  Figure  1. 


(«  (A ( xy  ) y  ) 

(A* 

(«  (A(xy  )y  ) 

(A(xy )y) 

(A* 

(<  (A(xy)y) 

(A* 

('  ( A( xy  )y  ) 

(Aj  (x  { A ( xy  )x  ) J — L ) ) 

(A*  (*  (A(zy  )x  )j — )))) 

(A*  (•  (A(xy  )t )_! — !_)))))) 

(A. 

(«  ( A ( xy  )y  ) 

(A* 

(«  (A(xy  )y ) 

(A(xy  )y  ) 

(Ax 

(«  (A(xy  )y  ) 

(Ax 

(x  (A(xy  )y ) 

(Ax  (x  (A(xy  )x  )j — !_)) 

(Ax 

(x(A(xy)y) 

(Ax 

(x  (A(xy  )x  )j — L ) ) 

(A. 

(.(A(ry)x)xx)))))) 
(Ax  (x  (A(zy)t)±±)))))) 

(A,  (x  (A(xy  )x  )j — !_))))))) 

igure  1.  Normal  Form  of  1/2 


6.1.6  Nonstandard  Objects: 

There  are  many  equivalent  representations  of  the  same  rational  number;  alternately,  these  may  be  con¬ 
sidered  different  (nonstandard)  rational  numbers  that  behave  the  same  under  the  arithmetic  operations. 
For  example,  +1/+2,  +2/+4,  -1/-2,  -2/-4  etc.  all  behave  the  same.  Similarly,  -+-0/+1,  +0/-1,  +0/+2, 
-0/-1,  etc.  all  behave  the  same.  “Infinite”  rational  numbers,  such  as  +1/+0,  -1/-0,  etc.  are  also  permitted 
by  the  representation,  and  can  be  considered  nonstandard  objects.  Operations  on  nonstandard  objects  are 
discussed  further  at  the  end  of  this  section. 


6.1. 7  Extract  Numerator:  numg  hd 


Comments:  Note  that 


nuntg(ffl/  r)  «*  hd(m/n)  =*  hd  <m,n>  «*■  m 
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6.1.8  Extract  Denominator:  deng  =»  hd  •  tl 
Comment*:  Note  that 

den  q  (m/n)  =s>  hd  “  tl  (m  /  n  )  =»  hd  •  tl  <  m  ,n  >  =»  n 

6.3  Arithmetic  Operation* 

6.2.1  Sum: 

sum<j  =»  ratg  :[sumz  :(fl  prod z  numg  deng  ,  Q  prod z  deng  numg), 
fl  prodz  deng  deng  i 

Comment*:  Alternately,  replacing  fl  with  a  dyadic  composition  we  have 

sumg  =*  ratg  :{sumz  :[prodz  :(numg  0 1st,  deng  •  2nd),  prodz  :(deng  *  1st,  numg  ”  2nd) j 
prodz:(deng  •  1st,  deng  “2nd)} 

This  definition  makes  use  of  the  equation 

m  /  n  -r  m  ' /  n  '  =  (mn  ’  +  nm ')  /  nn  ' 

To  see  this  note  that 

sumg  (m  /  n)  (m  ‘  /  n  ') 

=9  ratg  :[sumz :  (IT  prodz  numg  deng,fiprodz  deng  numg), 
fl  prodz  deng  deng)  (m  /  n)  (m  /  n  ) 

=s>  ratg  [sumz:(fl  prodz  numg  deng,  fl  prodz  deng  numg)  (m  /  n)  (m  /  n  )] 

[(fl  prodz  deng  deng  )  (rr»  /  r» )  (rr»  /  n  )[ 

Now  note  that 

(fl  prodz  deng  deng )  (m  /  n)  (m  '  /  n  ) 

=»  prodz  [deng  (m  /  n  )]  deng  (m '/  n  '))  =*  prodz  n  n  =  nn 


Also  note  that 


sum/  (n  prod ^  numg  deng  ,  fl  prod/  den g  numg)  (m/  n)  (m  /  n  ) 

=*  sum/  [0  prod/  numg  den*,  (m  /  n  )  (m  '/  n  ')|  [0  prod/  deng  num,  ("»/■)  (m  7  n 
=*  sum/  {prod/  [nums  (m/  n  )]  deng  (m  ' /  n  );} 

{prod/  [deng  (m  /  n)  !num<j  (m  7  *•  )i) 

=»  sum/  (prod/  m  n  )  (prod/  n  m  ) 

=  mn  '  +■ nm 

6.2.2  v46so/ute  V^a/ue.  absg  =»  ratg  •  (abs/  “numg,  absz  “deng) 

Comment*:  This  definition  is  based  on 

j  m  /  n  |  =  |  rn  I  /  |n| 

To  see  this  note 

absg  (m/  n)  =*  ratg  “  (abs/  •  numg  ,  abs/  •  deng  )  (m  /  n  ) 

=»  ratg  |(abs/  •  numg  )  (m  /  n  )]  i(absz  "  deng  )  (m  /  n  )] 

=>  ratg  (abs/  m )  (abs/  n ) 

6.2. S  Negation:  negg  =9  ratg  •  (neg/ 0  numg  ,  deng  ) 

Comment*:  This  is  based  on  -  (m  /  n  )  =  (-  m  )/  n  .  To  see  this  note 

negg  (m/n)  =>  ratg  "  (neg/ “  numg  ,  deng  )  (m  /  n  ) 

=>  ratg  ((neg/  "numg)  (m  /  n)j  [deng  (m  /  n  )] 

=»  ratg  (neg/  m  )  n 

6.2.4  Difference:  difg  =»  sumg  :( 1st,  negg  °  2nd) 

Comment*:  Note  that 

difg  r  *  =*  sumg. (1st, negg  "2nd)  r  *  =9  =>  sumg  r  (negg  «) 

6.2.5  Product:  prodg  =»  ratg  :(*  prod/  numg  ,*  prod/  deng) 

Comment*:  This  is  based  on  the  identity 


(m  /  n)  (m  '  /  n  ')  =  (mm  ')/  (nn  ') 


Note  that 


prodg  (m  /  n  )  (m  '/  n  ')  =a  ratg  :(#  prodz  numg  ,  *  prodz  deng  )  (m  /  n  )  (m  7  /t  ') 

=?>  ratg  [(*  prodz  numg  )  (m  /  n  )  (m  '/  n  ')1  (*  prodz  deng  )  (m  /  n  )  {m  '/  n  ')) 

Now  note  that 

(¥  prod z  numg  )  (m  /  n  )  (m  ' /  n  ')  =*  prodz  [numg  (m  /  n  )j  (numg  (m  ' /  n  ')] 

=>  prod  2  m  m 

Similarly, 

(♦  prodz  deng  )  (m  /  n  )  (m  ' /  n  ')  =»  prodz  n  n 

6.2.6  Reciprocal: 

recipg  =s>  ratg  0  (deng  ,numg  ) 

Comment*:  Note  that 

recipg  (m  /  n)  =s»  ratg  "  (deng  ,nmng  )  (m  /  n  ) 

=»  ratg  [deng  (m  /  n)]'[numg(m/  n)|  =»  ratg  n  m  <=>  n/m 

Notice  that  the  reciprocal  of  zero  is  defined  (has  a  normal  form)  since  recipg  0/  1  <=»  -*-1/  +0,  although 
this  is  a  nonstandard  object. 

6.2.7  Quotient,  quog  =*  prodg  :(lst,  recipg  °  2nd) 

Comments:  Note  that 

quog  r  «  =*  prodg  :(lst,  recipg  “2nd)  r  *  =»  prodg  r  (recipg  «) 

6.2.8  Floor:  floorg  =»  quoz  •  (numg  ,deng  ) 

Comments:  The  floor  of  a  rational  number  is  the  greatest  integer  less  than  that  rational  number.  Thus 
the  floor  of  m  /  n  is  the  integer  part  of  the  quotient  of  m  and  n  : 

floor  g  (m/n)  =>  quoz  •  (numg  ,deng  )  (m  /  n  ) 

=*  quoz  (numg  (m  /  n  )j  [deng  (m  /  n  )  =» 


quoz  m  n 


6  2.9  Ceiling:  ceil q  =*  floor,}  0  sum,}  (1/  2)1 


Comments:  The  ceiling  of  a  rational  number  is  the  least  integer  greater  than  or  equal  to  that  rational 
number.  To  compute  this  we  use  the  relation 

ceil^  r  =  floorg  (r  +  i) 

Note  that 

ceilg  r  =>  floor,}  "[sum,}  (1/  2)|  r  =»  floor<j  Isumg  (1/  2)  r] 

6.5  Predicates  and  Relations 

6.3.1  Nonnegative:  nonneg?}  =>  equiv  •  (neg?^  ■  numg  ,  neg?*  •  den$  ) 

Comments:  A  rational  number  is  nonnegative  if  both  its  numerator  and  denominator  have  the  same 
sign.  That  is  m  /  n  $  0  if  and  only  if 

equiv  (neg?^  m)  (neg?^  n)  <=>  true 

Note  that 

nonneg?}  ( mj  n )  =>  equiv  °  (neg?^  °  numg  ,  neg?^  “den^  )  (m  /  n  ) 

=*  equiv  [(neg?2  •  numa  )  (m  /  n  )]  |(neg?z  °  den,;  )  (m  /  n  )j 
=»  equiv  (neg?^  m)  (neg?z  n) 

6  3.2  Negative:  neg?<j  =»  not "  nonneg?g 

Comments:  A  rational  is  negative  if  it  is  not  nonnegative. 

6  3.8  Less-Than:  ltg  =»  neg?Q:difg 

Comments:  A  rational  r  is  less  than  a  rational  s  if  and  only  if  r  -  s  is  negative: 

It q  r  s  =*  neg?g.difg  r  s  =*>  neg?9(diffl  r  s) 

6.8.4  Greater- Than:  gt$  =*  conv  Itg 

Comments:  Note  that  gt<j  r  s  =»  conv  gt$  r  s  =>  gt q  s  r 
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6.3.5  Less  or  Equal:  leg  =>  not:gtg 

6.3.6  Greater  or  Equal:  geg  =*  not:ltg 

6.3.7  Equality:  eqg  =*■  eqj:(flprod^  numg  deng  ,  Cl  prod^  deng  num9) 

Comments :  Alternately,  we  can  replace  the  fis  with  dyadic  compositions: 

eqg  =*  eq^  :(prodj  :(num(j  “1st,  deng  “2nd), 
prod^:(deng  0  1st,  numg  °2nd)| 

Two  rationals  m  /  n  and  m  ' /  n  are  equal  if  and  only  if  mn  =  m  '  n  .  Note 
eq0  (m  /  n)  (m  '  /  n  ') 

=>  eqz  :(fl  prod^  numg  deng  ,  fl  prod/  deng  numg  )  (m  /  n  )  (m  '/  n  ') 

=>  eq^  ((fl  prod^  numg  deng  )  (m  /  n )  (m  ' /  n  ')]  j(fl  prod^  numg  deng )  (m  /  n  )  (m  */  n  ')] 

=*■  eq^  {prod^  [numg  (m/  n  )|  [deng  (m  '/  n  ')]}  {prodz  [deng  (m  /  n )]  [numg  (m  '/  n  ')]} 

=>  eq^  (prod^  m  n  ')  (prod^  n  m  ') 

( (ge<?  *  V  )  -  * 

6.3.8  Maximum:  maxg  =*  A(z»)1ei8ey 

0.4  Operations  on  Nonstandard  Rational  Numbers 

6.4.1  Infinite  and  Indefinite  Rational  Numbers: 

oo  =*  -i-l/  +0 
—  OO  =s»  -1/  4-0 
X  =5“  -*-0/  -rO 

Comments:  The  chosen  representation  for  rational  numbers  permits  the  nonstandard  rationals  1/0, 
-1/0  and  0/0.  We  sh  ^w  later  that  other  nonstandard  rationals,  such  as  2/0  and  -3/0,  are  equal  to  one  of 
the  three  defined  above. 

In  this  section  we  explore  the  extent  to  which  these  nonstandard  numbers  can  be  manipulated  con¬ 
sistently,  and  the  extent  to  which  they  have  the  expected  properties.  Although  one  goal  of  this  report  is 
to  show  the  development  of  real  analysis  without  recourse  to  actual  infinities,  there  is  no  problem  with  the 
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computable  manipulation  of  finite  symbols  representing  infinities. 


6.4.2  Predicate*: 

nonstandard^  =*  zero?2*deng 
standard?^  =*>  not  •  nonstandard^ 

infinite?g  =*  and  •  (nonstandard^  ,  not 0  *ero?z  •  numg  ) 
indefinite?,}  =*  and  •  (nonstandard^  ,  zero?^  "  numg  ) 

Comment*:  Thus,  infinite?  00  «=*  true,  nonstandard^  00  «=»  true,  etc. 

6.4.6  Standard  Equality: 

eqg  im  /  n  )  00  «=►  *qg  (m  /  n)  ~°°  <=>  zero ?w  « 
eqg  00-00  «=>  eqg  —0000  <=>  true 

Comment*:  These  results  follow  from  the  definition  of  ‘eqg  ’,  and  reflect  the  fact  that  that  definition  is 
only  appropriate  for  finite  rational  numbers.  In  the  following  subsection  we  extend  the  equality  operation 
to  accommodate  nonstandard  rationals. 


6.4.4  Nonstandard  Equality: 


eq  <3- 


and  (infinite?  7)  (infinite?  r) 
else  eqg  g  r 


equiv  (pos?z  "numg  q)  (pos?^  °  numg  r) 


TABLE  1.  Equality  Relation  on  Nonscandard  Rationals 


r 

00 

-00 

X 

7 

eqg  q  r 

false 

false 

true 

00 

false 

true 

false 

true 

-OO 

false 

false 

true 

true 

X 

true 

true 

true 

true 

Comment*:  This  operation  has  the  properties  shown  in  Table  1.  Notice  that  the  indefinite  quantity  x 
is  equal  to  all  quantities,  finite  or  infinite.  Therefore  this  relation  is  not  transitive,  since  0/1  equals  *,  and 
X  equals  1/1,  but  0/1  is  not  equal  to  1/1.  On  the  other  hand,  the  infinities  (oq  -oc)  are  unequal  to  all 
finite  quantities  and  are  unequal  to  each  other.  Henceforth,  we  use  '=’  as  an  abbreviation  for  ‘eqg+\ 
That  is,  'q  =  r  ’  means  ‘eqg  +  q  r\ 
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6.4-5  Nonnegative: 


nonneg?<j  oo  *-*  true 
nonneg?(j  -oo  <=»  false 
nonneg?<j  x  true 

Comment*:  These  properties  follow  directly  from  the  definitions: 

nonneg?<j  oo  =*  nonneg?g(  +  l/  +0) 

=»  equiv  (neg?2  -rl)  (neg?^  +0) 

=s>  equiv  false  false 
<=»  true 

Similarly, 

nonneg?g  -oo  =»  equiv  (neg?z  -1)  (negz  +0)  <=*  false 
nonneg?<j  x  =»  equiv  (neg?z  +0)  (neg?2  -0)  «=»  true 

Comment*:  As  expected,  oo  is  nonnegative  and  -oo  is  negative.  That  x>  the  indefinite  quantity,  is 
nonnegative  may  seem  incorrect.  This  situation  is  a  result  of  our  considering  all  zeros  to  be  nonnegative, 
including  -0.  Otherwise  we  could  define  a  “negative  indefinite”  quantity  =»  -0/  +0.  It  can  be  seen 
that  it  is  reasonable  to  consider  x  nonnegative,  since  x  is  equal  to  every  rational  number,  including  zero. 
Hence,  it  is  nonnegative,  in  the  sense  that  it  is  greater  than  or  equal  to  zero. 

6-4-6  Negative: 

neg?g  oo  =»  false 
neg?<j  -oo  =9  true 
neg?g  X  =*  false 

6-4-7  Nomtandard  Addition: 

/and  (eqg+  q  r)  (and  (infinite?  q)  (infinite?  r)j  -*  q 
sumg+  =*  A (?r)|el8e9um9  9  r 

Comment*:  The  properties  of  this  operation  are  shown  in  Table  2.  These  follow  from  the  definitions, 
for  example 
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TABLE  2.  Addition  of  Nonstandard  Rationals 


q  |  q  -r  r  X  -  X  X 

oo  j  oo  oo  X  X 

-oo  !  -oc  X  X 

X  I  X  XXX 

m/n+1/0  =  (mO  +  nl)/  |nO]  =  m/0=  oq  if  m  >  0 

m  /  n  +  -  1/0  -  (m  -0  +  n  ■— 1)/  (n  -0)  =  —  m/0  =  -  oq  if  m  >  0 

m  /  n  +0/0  =  (m  -0  +  n  -0)/  (n  0)  =  0/  0  =  jf 

A  nonstandard  addition  is  necessary  if  the  operation  is  to  have  the  expected  properties,  otherwise  we 
would  have  oo  +  oo  =  x,  since 

00+  00  =  1/  0  +  1/  0  =  (1  0  +  0  1)/  (00)  =  0/  0  =  x 

6-4-8  Negation: 

neg(j  oo  <=>  -oo 
negQ  -oo  <=>  oo 
"egg  X  <=*  X 

6-4  9  Difference: 

The  properties  of  subtraction  with  respect  to  the  nonstandard  quantities  follows  directly  from  those  of 
addition  by  the  identity 

difg  q  r  =  sum9  q  (negg  r) 

6-4-10  Ordering:  ltg+  =>  neg?g  :sum<j  +:(lst,negg  "2nd) 

TABLE  3.  Ordering  of  Nonstandard  Rationals 


lt(J  + 

r 

OO 

—  OO 

X 

q 

q  <  r 

true 

false 

false 

oo 

false 

false 

false 

false 

-oo 

true 

true 

false 

false 

false 

false 

false 

false 

false 

Comment $:  The  properties  of  this  operation  are  shown  in  Table  S.  As  expected  we  have  -oo  <  q  <  oo 
by  this  ordering.  Notice  that  x  '*  neither  greater  than  q  nor  less  than  q ,  for  it  is  equal  to  q ,  for  every  q  , 
finite  or  infinite. 
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prod0 

0 

r 

00 

-oo 

X 

0 

0 

0 

X 

X 

X 

<? 

0 

V 

ioo 

+  O0 

X 

oo 

X 

xoo 

00 

-oo 

X 

-oo 

X 

+  00 

—  oo 

00 

X 

X 

X 

X 

X 

X 

X 

The  properties  of  multiplication  of  nonstandard  rationals,  as  deduced  from  the  standard  definition  of 
rational  multiplication,  are  shown  in  Table  4.  The  sign  ±  must  be  taken  in  the  same  sense  as  the  sign  of 
q  (or  r  ),  whereas  the  sign  +  must  be  taken  in  the  opposite  sense. 

The  nonstandard  properties  of  the  reciprocal  are: 

recipe  oo  =  recipg  —  oo  =  0 

recipQ  x  =  X 

etc. 

The  properties  of  the  quotient  follow  from  the  identity 

quo9  q  r  =  prodg  q  (recipQ  r ) 
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7.  Real  Number* 


In  this  chapter  we  reach  the  culmination  of  our  development,  defining  real  numbers  as  convergent 
computable  rational  sequences  of  a  particular  kind.  Operations  on  real  numbers  are  defined  as  comput¬ 
able  operations  on  these  sequences. 

7.1  Primitive  Ideas 

7.1.1  Definition  of  a  Regular  Sequence: 

A  sequence  of  rational  numbers  z  is  called  regular  if  and  only  if  for  any  nonzero  natural  numbers  m 
and  n  we  have 

leg  [absg  (difg  zm  z„  )j  [suing  (  +  1/  m  )  (+1/  n  )j  «=*  true 
Intuitively,  the  sequence  z  is  regular  if  and  only  if 

I  z*,  -  I  ^  m'1  +  n"1 

for  all  natural  m  ,  n  >0.  Note  that  for  any  given  m  and  n  it  is  finitely  decidable  whether  this  condition 
holds.  However,  it  is  not  in  general  finitely  decidable  for  arbitrary  m  and  n  . 

7.1.2  Definition  of  Rea!  Number: 

A  real  number  is  a  regular  sequence  of  rational  numbers. 

7.1.3  Equality: 

Two  real  numbers  z  and  y  are  equal,  written  z  =  y ,  if  and  only  if  for  any  nonzero  natural  number  n 
we  have 

leg  [absg  (difg  z„  yB)|  [prodg  (+2/  n )]  <*==>  true 
Comment t:  Intuitively,  z  =  y  if  and  only  if 

I  *»  ~  V.  I  2»_1 

for  all  natural  n  >  0. 

Note  that  while  this  condition  is  finitely  decidable  for  a  given  n  ,  it  is  not  generally  finitely  decidable 
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for  arbitrary  n . 


7.1.4  Nonstandard  Objects:  There  are  two  alternative  notions  of  real  numbers.  First,  real  numbers  can 
be  considered  the  equivalence  classes  under  the  equality  relation  of  regular  sequences  of  rational  numbers. 
Second,  real  numbers  can  be  identified  with  the  regular  sequences  themselves.  In  this  case  we  find  that 
there  are  many  different  (nonstandard)  real  numbers  that  are  treated  as  equal  by  the  arithmetic  opera¬ 
tions  and  the  equality  relation. 

We  can  also  consider  the  nonstandard  real  number  all  of  whose  rational  approximations  are  oa  To  see 
if  this  object  Fits  the  definition  of  a  real  number  we  must  see  if  the  difference  between  its  m  and  n  th 
rational  approximations  is  greater  than  or  equal  to  1/  m  +  1/  n  .  Now,  the  difference  between  its  m  th 
and  n  th  approximations  is: 

!<*-<»!  =  I  x!  =  x 

We  have  previously  seen  that  x  is  greater  than  or  equal  to  every  rational  number,  since  it  is  equal  to 
every  rational  number.  Thus, 

|  oo-  oo]  =  x  ^  1/  m  +  1/  n 

and  so  the  sequence  composed  of  all  oa  is  a  real  number.  The  same  applies  for  -oa  The  following  opera¬ 
tions  work,  or  can  be  extended  to  work,  on  infinite  reals  in  much  the  same  way  as  the  rational  operations. 

7.1.5  Creation  of  Reals  from  Rationals:  reals  =*  Az(Anz) 

Comments:  This  function  converts  a  rational  number  into  a  real  all  of  whose  rational  approximations 
are  the  same: 

real*  r  =*■  Az(Anz)r  =>  Anr 
Now  notice  that  the  k  th  rational  approximation  of  real*  r  is: 

(real*  r)t  =»  real*  r  k  =*  Anr  k  =*■  r 

7.1.6  Real  Denotations: 
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0.0  =%  realR  (0/  1) 
1.0  =>  reals  (1/  1) 


7.1.7  Example  Normal  Form:  For  concreteness  we  exhibit  the  normal  form  of  the  real  number  0.0: 

0.0  =»  real*  (0/  0) 

=»  An  (0/0) 

=9  An(A«(j  A  ( xy )  y  A*  (<  A(zy)y  A(zy)y  A*  («  A(zy)y  As  («  A(zy)z  j_  j_) 

As  («  A ( zy  )  z  x  j_))) 

As  (*  A(zy)y  As  (*  A(zy)y  A(zy  )y  As  (s  A(zy)y  As  (s  A(zy)z  j_  j_) 

As  (s  A(zy)z  x  J_)))  As  (s  A(zy  )z  x  J_))) 

7.1.8  Example  Normal  Form  of  Irrational  Number:  We  have  said  that  this  computational  approach  to 
real  analysis  represents  all  objects  as  finite  structures.  To  illustrate  this  we  will  show  the  finite  structure 
that  represents  an  irrational  number,  the  square  root  of  2.  The  square  root  of  two  will  be  a  regular 
rational  sequence  s  such  that  st  is  the  Jbth  rational  approximation  of  the  square  root.  By  Newton’s 


Method  observe  that 


*4  -  -j(2  /  *4-1  +  *4-l)>  fc  >  1 


The  initial  approximation,  st  is  1.  Thus  the  sequence  can  be  defined  recursively  by  the  rule: 


rec  s  =  A 


k  -  1  -»  1 

el®e  *4-1  +  *4  -l) 


Since  this  sequence  converges  quadratically  it  is  regular  and,  hence,  represents  a  real  number.  To  show 
that  «  is  indeed  a  finite  formula  in  the  lambda  calculus,  we  rewrite  it  using  only  the  notational  abbrevia¬ 


tions  introduced  so  far: 


( (*q*  It)  -  1/  I 

*  “  [prodQ  (1/  2) 


(sum„  (quo„  (2/  1)  [t  (pred*  *)]}[»  (pred*  k] 


The  above  formula  is  finite;  it  will  remain  finite  when  all  the  abbreviations  are  eliminated,  since  each 


symbol  is  rewritten  by  a  finite  string.  We  illustrate  the  first  few  steps  of  this  process  (we  have  replaced 
the  name  ’  by  S/2’): 


( (eqjv  It)  -  1/1 

^  Y  X*  A*\else  [prodg  (1/  2)  (su mg  {quog  (2/  1)  j«  (predy  t)!}  (s  (predy  t)j)] 

=*  Y  As  At  {(eqy  1  t)  (1/  1)  [prod<j  (1/  2) 

(suing  {quog  (2/  1)  Is  (predy  *  )]}  [«  (predy  t )!)] } 

=»  Y  As  \k  {(not:ney  1  t)  (ratg  1  1)  [prodg  (ratg  1  2) 

(sumg  {quo<j  (ratg  2  1)  [a  (predy  k)\ }  [a  (predy  t)])]} 

By  eliminating  all  the  abbreviations  we  would  arrive  at  a  seminormal  form  expression  representing  the 
square  root  of  2.  The  full  seminormal  form  of  the  square  root  of  2  is  shown  in  the  Appendix. 

7.1.9  Rational  Approximation:  approx#  =>  conv  Id 

Comments:  The  function  ‘approx#  n  ’  gives  then  nth  rational  approximation  of  a  real  number.  That 
is,  if  x  is  a  real  number  then  ‘approx#  n  x  ’  is  its  n  th  rational  approximation: 

approx#  n  x  =*  conv  Id  nr  =»  Id  2  n  =»  z  n  «=*  xn 

7.3  Arithmetic 

7.2.1  Sum:  sum#  =>  A(zy)[sumg  “(z,y)  *  (prod*  2}j 

Comments:  This  definition  makes  the  nth  rational  approximation  of  z+y  the  sum  of  the  2nth 
rational  approximations  of  x  and  y  : 

(*  +  v)B  = 

Now  note  that: 

(sum#  x  y)„  =>  sum#  r  y  n 

=*  [sumg  "(z,y)  0  (prod*  2)]  n 
=*  sumg  *  (z  ,y )  ((prody  2)  n] 

=*  sumg  •  (z  ,y  )  (prody  2  n  ) 

=9  sumg  [z  (prody  2  n)|  [y  (prody  2  n  )j 

It  is  easy  to  show  that  z+y  is  a  real  number  and  that  this  operation  has  the  expected  properties  (commu¬ 
tativity,  associativity,  etc.). 
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7 .2.2  Maximum:  max* 


=>  $  maxg 

Comment*:  The  n  th  rational  approximation  of  the  maximum  of  two  reals  is  the  maximum  of  their 
n  th  rational  approximations: 

(max*  x  y  )n  =>  maxj  t  y  n  =*  $  maxg  x  y  n  =>  maxg  z„  yn 
7.2.8  Negation:  neg*  =>  B  negg 

Comment *:  The  n  th  rational  approximation  of  the  negation  of  a  real  is  the  negation  of  its  n  th 
rational  approximation: 

(neg*  x)n  =>  neg*  x  n  =>  B  negg  x  n  =*  negg  (z  n)  «=»•  negg  x„ 

7.2.4  Minimum:  min*  =»  neg*  :(♦  max*  neg*  ) 

Comment*:  This  is  based  on  the  fact  that  the  minimum  of  z  and  y  is  the  negation  of  the  maximum  of 
the  negations  of  x  and  y  : 

min*  z  y  =»  neg*  :(*  max*  neg*  )  x  y 

=»  neg*  {#  max*  neg*  z  y)  =>  neg*  (max*  (neg*  r)(  neg*  y)J 

7.2.5  Absolute  Value:  abs*  =»  S  max*  neg* 

Comment *:  The  absolute  value  of  a  number  is  the  maximum  of  the  number  and  its  negation: 
abs*  z  =»  S  max*  neg*  z  -*  max*  z  (neg*  z) 

7.2.6  Canonical  Bound:  cbd*  =»  ceilg  •  jsumg  (2/  1)J  "  absg  (approx*  1) 

Comments:  The  canonical  bound  of  a  real  is  used  in  the  definition  of  the  product,  which  follows.  The 
canonical  bound  of  z  is  the  least  integer  that  is  two  greater  than  the  absolute  value  of  the  first  rational 
approximation  of  x  : 

cbd*  x  =*  {ceilg  •  Jsumg  (2/  1)]  *  absg  (approx*  1)}  x 
=*  ceilg  {sumg  (2/  1)  [absg  (approx*  1  x)]} 

=*  ceilg  jsumg  (2/  1)  (absg  x ,) j 

The  latter  represents  the  ceiling  of  2  +  j  xt  |  . 
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7. 2.7  Product: 


prods  =9  A(xy)|prod,}  ‘(i,k)  °  (prod/y  2)  •  mag/  •  .prod/  (♦  max/  ebd*  x  y)j| 

Comments:  The  n  th  rational  approximation  of  xy  is  the  product  of  the  2 kn  th  rational  approximations 
of  z  and  y  ,  in  which  k  is  the  maximum  of  the  canonical  bounds  of  x  and  y  .  First  notice  that 

®  max/  cbdR  x  y  =*  max/  (cbdg  x)  (cbd*  y)  <=*■  k 

Therefore, 

(prod/}  x  y)„  =»  prods  x  y  n 
=»  [prod(j  °(x,y)  •  (prod/y  2)  •  mag/  •  (prod/  i)j  n 
=*  prod,}  “  (x ,y )  {prod/y  2  [mag/  (prod/  k  n  )]} 

=*>  prod}}  •  (x  ,y  )  (2tn  ) 

=s>  prodg  (x  2 kn)  (y  2kn  ) 

<f=»  prodg  X  2kn  V  2tn 

7 .2.8  Exponentiation  to  Natural  Power:  expts  =»  A(nx  )[rpt  (prod/}  x)n  1] 

Comments:  The  basis  for  this  is: 

expt8  n  x  =*  rpt  (prod/}  z)  n  1  =9  (prod/}  x )"  1  =  x" 

We  usually  write ‘x"  ’  for ‘expt/}  n  x 

7 .2.9  Notation : 

Since  each  of  the  classes  of  objects  we  have  defined  effectively  includes  the  preceding,  and  the  opera¬ 
tions  and  relations  of  the  latter  classes  extend  the  former  classes,  we  will  drop  the  subscripts  (/V,  Z ,  Q 
and  R  )  and  functional  notation,  and  henceforth  use  conventional  mathematical  notation.  Thus  we  will 
write  ‘x  •+•  2y  ’  for 

sum/}  x  {prod*  (real*  (2/  l)j  y } 
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7.S  Predicates  and  Relations 

7.5.1  Positive  with  Modulus:  pos?  =*  A(nx)!gt<3  xH  (1/  n)l 

Comments:  If  n  is  a  nonzero  natural  number  and  z  is  a  real  number,  then  z  is  called  positive  with 
modulus  n  if  and  only  if  xn  >1/  n  .  The  argument  n  names  a  place  where  zn  deviates  sufficiently  from 
zero.  Notice  that  this  predicate  is  algorithmically  decidable;  that  is,  if  z  and  n  are  defined  then 
pos?  n  x  is  defined. 

7.8.2  Positive:  We  say  that  a  real  number  is  positive  if  and  only  if  there  is  a  nonzero  natural  n  such  that 
the  real  is  positive  with  modulus  n  .  That  is,  x  is  positive  if  and  only  if  there  is  a  nonzero  natural  n  such 
that  pos?  n  i  <=»  true. 

Comments:  Note  that  it  is  not  in  general  algorithmically  decidable  whether  or  not  a  real  number  is 
positive.  Thus,  it  is  not  safe  to  treat  a  real  number  as  a  positive  unless  we  have  calculated  (i.e.,  con¬ 
structed  algorithmically)  the  modulus  n  required  by  the  above  definition.  This  algorithmic  undecidability 
is  inherited  by  all  the  other  relations. 

7.S.S  Negative  with  Modulus:  neg?  =s>  pos?:(lst,  neg  *  2nd) 

Comments:  Here  z  is  negative  with  modulus  n  if  and  only  if  -z  is  positive  with  modulus  n  : 
neg’  n  z  =*  pos?:(lst,  neg  •  2nd)  n  z  =>  pos?  n  (neg  x  ) 

Notice  that  if  -z  is  positive  (with  modulus  n  ),  then  (-z)„  >  1/  n  .  Therefore,  i  is  negative  (with 
modulus  n  )  if  and  only  if  z„  <  -  1/  n  . 

7 .8.4  Zero  with  Modulus  Predicate:  zero?  =»  not  •  or  *  (pos?, neg7) 

Comments:  A  number  is  zero  (with  modulus  n  )  if  and  only  if  it  is  neither  positive  nor  negative  (with 
modulus  n  ): 

zero'5  z  =»  not  •  or  •  (pos?, neg?)  z 
=>  not|or  •  (pos?, neg?)  z  j 
=*  not  or  (pos9  z  )  (neg?  z  )[ 

Notice  that  z  is  zero  with  modulus  n  if  and  only  if-l/i»$xII5  1,n,  that  is,  z„  $  1/  r» . 
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7 .3.5  Equality  with  Modulus:  eq  =>  An ’(zero?  n  ):dif; 

Comments:  Two  reals  are  equal  with  modulus  n  if  and  only  if  their  difference  is  zero  with  modulus  n  : 
eq  n  z  y  =»  (zero?  n  ):dif  z  y  =*  zero?(dif  z  y) 

Notice  that  since 

I  *  ~  V  I  n  =  !(*-»)«!  =  z2n  -  l/2» 

we  know  that  x  equals  y  with  modulus  n  if  and  only  if  j  z2n  -  yln  |  ^  1/  n  .  We  will  generally  write 
‘z  =„  y’  to  mean  that  z  equals  y  with  modulus  n  .  Finally,  notice  that  z  equals  y  if  and  only  if  for  all 
nonzero  naturals  n  ,  z  =  „  y  . 

7 .3.6  Sonnegative:  A  real  number  z  is  called  nonnegative  if  and  only  if  for  all  nonzero  natural  numbers 
n  we  have  z„  ^  -  1  /  n  . 

Comments:  This  is  also  in  general  algorithmically  undecidable,  since  it  depends  on  an  infinite  process, 
viz.,  checking  that  the  above  condition  applies  for  all  n  . 

7.3.7  Greater- Than  with  Modulus:  gt  =*  An  [(pos?  n  ):dif] 

Comments:  We  say  that  z  is  greater  then  y  with  modulus  n  when  z  -  y  is  positive  with  modulus  n  : 
gt  n  z  y  --=»  (pos?  n  ):dif  z  y  =*  pos?  n  (dif  z  y ) 

Hence,  since 

(*  -  v)»  =  z2n  -  V2n 

we  have  z  is  greater  than  y  with  modulus  n  if  and  only  if  z2n—  yjB>l/n,  that  is, 
Zj „  >  y2n  +  1/  n.  Notice  that  this  relation  is  algorithmically  decidable,  but  that  it  requires  supplying 
the  modulus  n  ,  which  names  the  place  where  z  sufficiently  exceeds  y .  We  will  generally  write  z  >„  y 
to  mean  that  z  is  greater  than  y  with  modulus  n  , 

7 .3.8  Greater-  Than  Relation:  We  say  z  is  greater  than  y  and  write  z  >  y  if  and  only  if  z  -  y  is  posi¬ 
tive. 

Comments:  That  is,  if  and  only  if  there  is  a  nonzero  natural  number  n  such  that  z  is  greater  than  y 
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with  modulus  n  : 


gt  n  x  y  -==>  true 

This  relation  is  not  algorithmically  decidable. 

7.3.9  Leu*  Than  with  Modulus:  It  An  conv  (gt  n  )] 

Comments:  Observe 

It  n  z  y  =*  conv  (gt  n  )  z  y  =>  gt  n  y  x . 

Thus.  i2„  •+•  1/  n  <  y2„ 

7 .3.10  Less- Than  Relation:  We  say  z  is  less  than  y  and  write  z  <  y  if  and  only  if  y  >  z 

7.3.11  Greater- or-  Equal  with  Modulus:  ge  =s>  An  [or:(gt  n  ,  eq  n  )] 

Comments:  To  see  this,  observe: 

ge  n  z  y  =s>  or:(gt  n  ,  eq  n  )  z  y 

=>  or  (gt  n  z  y  )  (eq  n  z  y  ) 

Comments:  Note  that  i  is  greater  than  or  equal  to  y  with  modulus  n  if  and  only  if 
Zln  -r  1/  n  >  yJn. 

7.3.12  Greater- or-  Equal  Relation:  We  say  z  is  greater  than  or  equal  to  y  and  write  z  ^  y  if  and  only  if 
z  -  y  is  nonnegative. 

7.3.13  Less-or- Equal  with  Modulus:  le  =»  An  (conv  (ge  n  )J 

Comments:  Note  that  z  is  less  than  or  equal  to  y  with  modulus  n  if  and  only  if  z7„  $  y7„  +  V  »  • 

7 .3.14  Less- or- Equal  Relation:  We  say  z  is  less  than  or  equal  to  y  and  write  z  $  y  if  and  only  if 
z  ^  y 

7.3.15  Inequality  with  Modulus:  ne  An  jnot:(eq  n  )j 
Comments:  Observe  that 

ne  n  z  y  =*  not:(eq  n)  z  y  =»  not(eq  n  z  y) 
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Comment s:  Note  that  x  is  not  equal  to  y  with  modulus  n  if  and  only  if  j  x2n  -  y  2n  j  >  1/  n  . 

7 .3.16  Inequality  Relation:  We  say  z  is  unequal  to  y  and  write  z  #  y  if  and  only  if  either  z  >  y  or 
z  <  y. 


7.4  Reciprocals  and  Quotients 

7.4.1  Reciprocal:  If  z  and  y  are  real  numbers,  then  z  is  called  a  reciprocal  of  y  if  and  only  if  zy  =  1. 

Comments:  The  following  three  sections  define  a  function  for  computing  the  reciprocal  of  a  real 
number. 

7.4.2  Zero  Bound:  zerobnd  =s>  A(nx)[floorc  (recip<j  {quo<}  [dif<j  (absg  z„)  (1/  n  )j  (2/  1)})] 

Comments:  If  z  is  nonzero  with  modulus  n  ,  then  zerobnd  n  z  has  the  following  property.  Let 
TV  =  zerobnd  n  z.  Then,  for  all  m  ^  TV,  zm  ^  1/  TV. 

7.4.8  Reciprocal  with  Modulus: 


recip 


( k  <TV 


A(nz  )  A k' 


k>N 


1/  (*„») 
V  I'm') 


where  TV  =  zerobnd  n  z 


Comments:  Notice  that  the  reciprocal  is  computable  only  if  a  modulus  of  nonnullity  is  known.  Intui¬ 
tively, 

(recip  n  z  )*  =  1/  (x^,),  if  k<N 

=  1  /(**,),  if  *2* TV 

We  will  generally  write  ‘z-1’  for  ‘recip  n  z’,  omitting  the  modulus  of  nonnullity.  It  is  important  to 
observe  that  z'1  is  computable  only  if  we  can  compute  this  modulus. 

7.4  4  Quotient  with  Modulus:  quo  An  {prod:[l«t ,  (recip  n  ) 0  2nd!} 

Comments:  Therefore,  the  quotient  is  the  product  of  the  first  and  the  reciprocal  of  the  second: 
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quo  x  y  =9  prod:'lst,  (recip  n  )  0  2nd;  x  y  =*  prod  x  (recip  n  y  )  =  xy-1 
We  will  generally  write  ‘x  /  y  ’  or  ’  for  ‘quo  n  x  y  Notice  that  x  /  y  is  computable  only  if  a  modulus 

y 

of  nonnullity  for  y  is  known.  The  modulus  is  implicit  in  the  division  notation. 

7.5  Important  Properties  of  the  Real  Numbers 

7.5.1  Approximation  Theorem:  The  n  th  rational  approximation  of  a  real  x  approximates  x  to  within 
1/  n  .  That  is,  for  all  nonzero  naturals  n  , 

I  *  “  *•  I  ^  1/  » 

7.5.2  Rational  Between  Two  Reals:  meanrat  =»  A(nxy)jprodg  -j  (suing  xlB  yjB)] 

Comments:  Notice  that  if  x  and  y  are  reals  and  n  is  a  nonzero  natural  number,  then  the  result  of 
(meanrat  n  x  y  )  is  a  rational  number  with  the  value 

meanrat  n  x  y  <=>  -|(x2b  +  y2„) 

Since  (meanrat  n  x  y )  is  the  average  of  the  2nth  rational  approximations  of  x  and  y ,  it  is  between 
these  approximations.  We  will  generally  write  ‘meanrat,  (x  ,y  )’  for  ‘meanrat  n  x  y’. 

7.5.5  Denseness  of  Rational  Numbers:  If  x  and  y  are  real  numbers  with  x  <  y,  then  there  exists  a 
rational  number  r  such  that  x  <  r  <  y . 

Proof:  If  x  <  y  then  there  is  some  nonzero  natural  n  such  that  x  <  y  with  modulus  n  .  We  now 
claim  that 

r  =  meanrat,  (x  ,y  ) 

is  a  rational  number  satisfying  x  <  r  <  y .  That  r  is  a  rational  follows  from  the  definition  of  meanrat. 
To  show  x  <  r  we  must  show  r  -  z  is  positive.  Therefore,  recalling  that  (x  -  y  )„  =  x2„  -  yJo,  we 
have 

r-x  ^  r  -  x2n 

r  ~  -  I  *2«  ~  *  I 

=  +  V2»  )  -  *2»  -  I  *2,  -  *  ! 


-47- 


=  \{Vin  ~  *2»  )  -  !  *J«  ~  *  | 

By  the  approximation  theorem: 

r  —  z  ^  "2(1^211  —  *2n)  —  (2n  ) 

Since  z  <  y  with  modulus  n  ,  we  know  y2n  -  zJn  >n~1,  therefore 

r  -  1  >  in'1-  (2n  )_1  =  0 

Therefore  z  <  r  .  A  similar  derivation  shows  that  r  <  y  . 

7.5.4  Cantor’s  Theorem:  Let  R  be  any  sequence  of  rational  numbers  and  u  and  t  any  real  numbers 
with  u  <  v  .  Then  there  exists  a  real  number  (  between  z  and  y,  u  <  (  ^  v ,  such  that  for  all  nonzero 
naturals  n  we  have  £  /  Rn . 

Proof:  Since  u  <  c  there  is  some  nonzero  natural  number  k  such  that  is  the  modulus  of  this  relation, 
u  <K  v .  We  will  construct  sequences  of  rationals  z  and  y  such  that  the  following  properties  hold: 

1.  For  all  natural  m  ^  n  $1, 

»  ^  zn  ^  zm  <  ym  ^  yn  ^  v 

That  is,  the  zn  s  are  less  than  the  yn  s,  and  they  are  always  growing  closer  together. 

2.  For  all  nonzero  natural  n  ,  either  z„  >K  Rn  or  yn  <.K  Rn  (or  both).  That  is,  the  x„  s  and  y„  s  are  all 
sufficiently  different  from  the  corresponding  Rn  s. 

3.  For  all  nonzero  natural  n,  yn  —  z„  <  1/  n.  That  is,  z  and  y  are  equal  (in  the  sense  of  real 
numbers).  We  will  define  £  =  z  =  y . 

W'e  construct  the  sequences  z  and  y  using  a  “diagonalization”  function  ‘diag’,  which  constructs  the 
two  sequences  in  parallel.  If  5  =  diag(*,u,c),  then  z  ~  S0  and  y  =  S,.  Hence  we  introduce  ‘x’  and  ‘y’ 
as  explicit  abbreviations: 

‘x  ’  =»  ‘diag  k  u  v  0’ 

‘y  ’  =*  ‘diag  k  u  v  1’ 

The  definition  of  the  diagonalization  function  is: 
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rec  diag  (k,u  ,v  ,b  )  =  An 
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$  max(rt„,  r„_,,  y  -  1/  n  )  <  r„  <  yB 

(Notice  that  in  this  case  x„  is  constrained  to  be  within  1/  n  below  y„  .)  Combining  these  results  we  have 

u  =  z0  <  z {  <  <*»<»*<  <  Vi  <  Ito  =  v 

Thus  condition  (1)  is  satisfied.  Also  notice  that  since  xB  and  yB  are  always  constrained  to  be  within  1/  n 
of  each  other,  y„  -  x„  <1/  n  ,  and  condition  (3)  is  satisfied. 

To  see  that  (2)  holds,  consider  any  step  in  the  construction  of  z  and  y.  Since  zn.x  <  y„_,,  either 
Rn  >„  x„-,  or  Rn  <„  yB_|  (or  both).  If  Rn  >K  xB_t  then 

z„  <  min(JiB ,  yB_,)  ^  Rn 

so  x„  <  Rn  ■  Similarly  yB  <Rn  On  the  other  hand,  if  J?B  <B  y„_i  then 

Rn  $  max(fl„,  xB_,)  <  yB 

so  Rn  <  yB .  Similarly  R„  <  xB  .  Thus  we  have  that  either  xB  >  Rn  or  yB  <  Rn . 

The  rest  of  the  proof  is  simple.  First  we  must  show  that  x  and  y  are  real  numbers.  Hence,  suppose 
that  m  $  n  .  By  properties  (I)  and  (S)  we  have 

I  “  *n  I  =  xm  ~  <»»"*«  <  1/  «  <  1/  m  +  1/  fl 

Thus  x  is  a  real  number;  the  result  is  analogous  for  y.  Furthermore,  by  (S)  we  know  that  x  =  y. 
Therefore,  let  £  =  x  =  y  . 

To  complete  the  proof  we  must  show  that  the  diagonal  number  {  does  not  appear  in  the  sequence  R  . 
That  is,  for  all  nonzero  natural  n  ,  £  #  R„  .  Let  n  be  chosen.  By  (2)  either  x„  >  Rn  or  y„  <  Rn .  Sup¬ 
pose  z„  >Rn  Then,  since  (1)  implies  x  >  x„  we  know  x  >  Rn  and  hence  {  =  z  ?  Rn .  On  the  other 
hand  suppose  that  y„  <  Rn  .  Since  (1)  implies  y  <  y„  we  know  y  <  Rn  and  hence  £  =  y  #  /fB .  So,  in 
either  case  £  #  Rn  ,  and  the  theorem  is  proved. 

Comments.  Cantor’s  Theorem  is  usually  interpreted  as  saying  that  the  real  numbers  cannot  be 
enumerated.  Indeed,  we  have  shown  here  how  given  any  enumeration  R  of  (computable)  real  numbers 
we  can  compute  a  real  £  not  in  that  enumeration.  But  this  is  paradoxical,  since  we  clearly  can  enumerate 
all  lambda  calculus  formulas:  just  enumerate  all  finite  strings  of  symbols  in  the  (finite)  alphabet,  and 
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strike  out  those  that  do  not  fit  the  syntactic  rules  of  the  lambda  calculus.  Since  this  enumeration  will 
include  all  those  formulas  that  compute  real  numbers,  it  seems  that  we  have  the  enumeration  of  the  (com¬ 
putable)  reals  that  Cantor’s  Theorem  denies.  This  is  in  effect  the  famous  Lbwenheim-Skolem  Paradox. 

We  can  escape  from  the  paradox  by  looking  carefully  at  the  premises  of  Cantor’s  Theorem,  which 
require  that  we  be  given  a  sequence  R  of  real  numbers.  Recall  that  a  sequence  of  real  numbers  is  defined 
as  a  lambda  calculus  formula  R  such  that  for  every  nontero  natural  number  n  ,  Rn  is  defined  (i.e.,  has  a 
seminormal  form),  and  Rn  is  a  real  number.  Now,  our  enumeration  of  lambda  calculus  formulas  is  not  a 
sequence  of  reals,  since  it  generates  many  formulas  that  do  not  fit  the  definition  of  a  real  number,  and 
may  not  even  have  a  seminormal  form.  Furthermore,  we  cannot  convert  our  enumeration  of  lambda  cal¬ 
culus  formulas  into  a  sequence  of  reals  by  striking  out  those  formulas  that  are  not  reals,  since  this  pro¬ 
perty  is  algorithmically  undecidable. 

This  leads  us  to  the  following  interpretation  of  Cantor’s  Theorem:  It  is  impossible  to  define  a  comput¬ 
able  enumeration  that  generates  all  the  computable  reals  and  only  the  computable  reals. 


% 
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8.  Calculus 


8.1  Convergence,  Limit*  and  Serie* 

8.1.1  Convergence: 

If  X  is  a  sequerce  of  reals  and  N  is  a  sequence  of  natural  numbers,  then  we  say  that  X  converge*  to  y 
with  modulus  N  if  and  only  if  for  all  nonzero  natural  numbers  k  and  n  such  that  n  ^  Nt  we  have 

I  x»  -  V  I  < 

When  N  or  y  is  not  relevant  we  may  omit  mentioning  them  and  say,  for  example,  X  converge*  or  X  i* 
convergent. 

8.1.2  Limit: 

limy  X  =>  S  [X  •  (limmod  JV)]  (prod  2) 
where  limmod  =>  A N  {max  •  [prod  3,  N  ‘  (prod  3)| 

Comment*:  Intuitively,  we  have: 

limy  X  =  At  l(ATMt)n  !  where  Af  =  At  |max(3i ,  Nik )] 

That  is,  the  k  th  rational  approximation  of  limy  X  is  the  2Jbth  rational  approximation  of  xMt  ,  where  the 

k  th  rational  approximation  of  M  is  the  maximum  of  Zk  and  Nu  .  To  see  that  our  definition  gives  this, 
let  M  =  limmod  N  and  observe: 

(limy  X  )t  =>  limy  Af  k 
=*  S  [Af  •  (limmod  JV ))  (prod  2)  k 
=»  \X  •  ( limmod  N)J  k  (prod  2  k) 

=?  X  (limmod  N  k  )  (prod  2  k ) 

=  X  [M  k)  (2*) 

=  {XMkhk 

Now  we  check  the  definition  of  ‘limmod’: 


limmod  N  k 


Mk  =* 

=»  max  •  [prod  5,  N  °  (prod  3)]  k 
=»  max  (prod  3  k )  \N  (prod  3  k  )j 
=  max  3  k  V  si 

Notice  that  the  limit  of  a  convergent  sequence  is  always  computable,  provided  that  the  modulus  of  con¬ 
vergence  is  supplied. 

8,1.3  Convergence  to  Limit* :  If  the  sequence  of  reals  X  converges  with  modulus  N ,  then  X  converges  to 
limy  X  (with  modulus  JV). 

Proof:  There  are  two  parts  to  the  proof.  First  we  must  show  that  limy  X  is  a  real  number;  second  we 
must  show  that  X  converges  to  limy  X  . 

Let  y  =  limy  X  .  To  show  that  y  is  a  real  number  we  have  to  show  that  for  any  nonzero  rationals  m , 
n  we  have 

i  Vm  ~  Vn  i  TO"1  +  *»"' 

Without  loss  of  generality  assume  that  m  ^  n  and  recall  that  M*  =  max(3 k,Ntk).  By  the  triangle 

inequality,  the  definition  of  y  and  the  approximation  theorem  we  have: 

:  Vm  -  Vn  i  ^  !  Vm  ~  XM^  |  +  |  XM^  -  XMm  |  +  |  XM*  -  yn  j 

=  I  ~  %Mm  I  +  I  !  +  i  fm  ~  I 

^  (2m)'1  -t-  )XUm-XU'\  +(2n)_1 

Since  X  converges  with  modulus  JV,  it  converges  to  some  real  number  with  this  modulus.  Call  this 
number  z .  Since  X  converges  to  z  with  modulus  N  we  know  that  for  any  n  ^  Nk  we  have 

'  X„  -  z  \  %  k~x.  Since  Mm  ^  N Sm  (by  the  construction  of  M),  we  know 

I  XUm  -  z  |  ^  (3m)'1 

Similarly  for  Mn  ^  N in  : 

I  XM%  -  z  |  $  (3n)-‘ 

Combining  these  two  inequalities  with  the  previous  inequality  yields 
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!  »m  -V*  !  ^  (2m  )  *  -  (3m)'1  -t-  (Sn  )  1  +  (2n)'‘ 

^  m  _1  +  n  _1 

Therefore  y  is  a  real  number. 

Now  we  must  show  that  .V  converges  to  y ;  in  fact  we  will  show  that  it  converges  with  modulus  M. 
Therefore,  suppose  that  n  ^  Af*  ;  we  must  show  that  y  -  j  <  k'1.  By  the  triangle  inequality,  the 
approximation  theorem,  the  definition  of  y  and  the  convergence  of  X  we  have: 

i  y  —  Xn  |  ^  I  y  -  yn  i  -1-  I  Vn  -  %m%  !  +  I  Xu 4  \ 

<  n  ' 1  +  |  ( -  XM *  |  +  |  XM -  Xn  | 

$  n'1  +  (2n  )_1  +  |  Xut~Xn  | 

$  «-*  +  (2 »)-«  *  !  Xm%-2  (  +  (  Xn-z  | 

Now,  since  Mn  >  N in  and  n  ^  M*  5  N St  we  have 

|y-*«|  ^  n-‘  +  (2n)-‘+  (3n)-‘ +  (3*)-1 

Now,  since  n  $  A/»  ^  3k  we  k now 

\y-Xn  I  s:  (3k)“l  +  (6k)'1  +  (9k)'1  +  (3k)'1 
=  (1/3  +  1/6+1/9+1/  3)k 'l 
=  (17/  18)k~ 1 

<  k'1 

Thus  X  converges  to  y  =  limy  X  with  modulus  M.  But  since  we  were  given  that  X  converges  to  its 
limit  with  modulus  N ,  we  know  that  X  converges  to  lim^  X  with  modulus  N . 

8.1-4  Propertiei  of  Lirrriti:  Let  X  and  Y  be  sequences,  then 

lim  sum  °  (X  ,  K)  =  sum(lim  X  ,  lim  K) 
lim  dif  *  (X ,  Y )  =  dif(lim  X  ,  lim  Y ) 
lim  prod  •  (X  ,  Y )  =  prod(lim  X  ,  lim  Y) 
lim  quo  "  ( X  ,  Y)  =  quo(lim  X  ,  lim  Y ) 

That  is,  the  limit  of  the  sum  is  the  sum  of  the  limits,  etc. 
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8.1.5  Compact  Interval:  A  pair  of  reals  'a  ,b]  is  called  a  compact  interval  if  and  only  if  a  $6  . 

Commenti:  When  we  say  that  a  real  z  is  in  the  compact  interval  a  ,6  j  we  mean  that  a  $  z  $6. 
Notice  that  it  is  not  in  general  algorithmically  decidable  whether  a  pair  of  reals  is  a  compact  interval,  or 
whether  a  real  is  in  a  given  compact  interval.  We  will  occasionally  use  other  interval  notations  with  their 
obvious  meanings.  For  example,  we  say  that  z  is  in  the  half-open  interval  (a  ,6  j  if  and  only  if 
a  <  x  ^  6  . 

8.1.6  Compact  Interval  with  Modulus:  A  pair  of  real  numbers  j  a  ,  6  ]  is  called  a  compact  interval  with 
modulus  n  if  and  only  if  a  is  less  than  or  equal  to  6  with  modulus  n  .  We  say  that  a  real  z  is  in  the  com¬ 
pact  interval  provided  that  a  is  less  than  or  equal  to  z  with  modulus  n ,  and  z  is  less  than  or  equal  to  6 
with  modulus  n . 

Comments:  If  the  modulus  n  is  provided,  then  it  becomes  algorithmically  decidable  whether  a  pair  of 
reals  is  a  compact  interval,  and  whether  a  given  real  belongs  to  the  interval. 

8.1.7  Notation  for  Limits:  When  the  modulus  of  convergence  jV  is  understood  or  implied,  we  will  write 
‘lim  X  ’  instead  of  ‘limN  X  ’.  We  also  permit  the  following  abbreviations: 

lim  /  (n  )  =»  lim  / 

r»  -»00 

lim  g  (n  )  lim  g  (-n  ) 

n  -*-00  »  -*00 

The  first  of  these  limits  exists  provided  that  there  are  N  and  y  such  that  for  all  k  >  0  and  all 

m  2  N(k)  we  have  /  (m  )  -  y  |  $  1/  t.  The  second  exists  provided  that  there  are  N  and  y  such 

that  for  all  Jt  <0  and  all  m  ^  N  (k)  we  have  |g(m)~yj  $  —  1/  k  . 

If  /  is  defined  on  (0,1]  and  g  is  defined  on  [-1,0)  then  we  can  write: 

lim  /  (z)  =s*  lim  /  (1/  n) 

x  -»0  t  n  -*oc 

lim  g{z)  =?«  lim  p(-z) 

« -o-  i-0  + 

The  first  of  these  exists  provided  that  there  are  6  and  y  such  that  for  all  e  >0  and  all  0  ^  f  ^  6  (« )  we 

have  /  (f )  -  y  \  $  e  .  The  existence  condition  for  the  second  is  analogous. 
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If  lim  /  (z  )  =  lim  /  { z ) ,  then  we  can  write  lim  /  ( * )  to  denote  their  common  limit. 

*  -0+  t  -*Q -  s  -*0 

If/  is  defined  on  (y  ,y  -*-l]  and  g  is  defined  on  [y  —  l,y  )  then  we  can  write: 
lim  /  (z  )  =*  lim  /  (a  -+-  6 ) 

t  —a  6-0-*- 

lim  y  ( z )  =»  lim  /  (a  6) 

t  -  6—0- 

The  first  of  these  exists  provided  that  there  are  6  and  y  such  that  for  all  c  >0  and  all  0  $  f  ^  6  (t)  we 
have  ;  /  (a  +  f)  —  y  j  $  c.  The  existence  condition  for  the  second  is  analogous.  We  use  lim  /  ( z  )  to 

%  — a 

denote  the  common  limit  when  the  limits  approaching  from  above  and  below  a  are  equal. 

If/  is  defined  on  the  compact  interval  I  =  [a  ,v  {  then  we  can  write: 

/lim  /  (z  )  =*  lim  /  ja  -t-6(v-a)] 

t  —a  *  4-°  + 

/  lim  /  (z )  =*  lim  /  [a  +  6  (a  — u  )) 

*  — a  - 

We  will  omit  mention  of  the  interval  /  and  the  direction  from  which  the  limit  is  approached  when  these 
are  not  relevant.  Again,  we  use  /  lim  /  (z )  for  the  common  limit. 

t  -  O 

/n=  0-0 

8.1.8  Partial  Sumt  of  Strict:  rec  E  X  =  An  i  else  S  X  (n  - 1) 

Comments:  A  series  is  a  sequence  that  is  meant  to  be  summed.  If  X  is  a  series  then  EX  is  the 
sequence  of  partial  sums  of  X  .  That  is: 

(EX)0  =  0 
(EX),  =  X, 

(EX),  =  X,  +  X, 

(EX).  =  X„  +  X„_,  +  •  •  •  -r  X,+  X, 

06 

8.1.9  Infinite  Sums  of  Series:  =»  lim  •  E 

00 

Comments:  If  X  is  a  convergent  series,  then  X  is  the  sum  of  all  its  terms. 
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Y.X  =>  lim-EX  =>  lim(EX)  =  Hm  (EX). 

n  —oo 

The  modulus  of  convergence  here  is  implicit;  if  we  must  be  explicit  about  it  we  will  write: 

0C 

£*  =*  lim*  •  £ 

00 

8.1.10  Exponential  Function:  exp  =»  Az  quo  "  (conv  expt  z  ,  fac)] 

00 

Comment*:  The  basis  for  this  definition  is  exp  x  =  An  (z*  /  n  !) 

8.2  Continuous  Functions 

8.2.1  Continuity  on  a  Compact  Interval:  If  /  and  u  are  real-valued  functions  and  /  is  defined  on  a 
compact  interval  / ,  then  we  call  /  continuous  on  /  with  modulu*  of  continuity  w  if  and  only  if  for  each 
positive  «  we  have 

!/(*)-  /  (if) I 

whenever  z  and  y  are  in  /  and  z  -  y  |  $  u;(e  ). 

Comment*.  We  often  omit  mention  of  the  modulus  of  continuity  when  it  is  not  relevant.  Notice  that 
it  is  not  in  general  algorithmically  decidable  whether  a  function  is  continuous,  since  that  decision  would 
require  testing  the  above  condition  for  every  r  and  y  in  the  interval. 

8.2.2  Continuity  on  an  Arbitrary  Interval:  We  say  that  a  real-valued  function  defined  on  an  arbitrary 
interval  is  continuous  on  that  interval  provided  that  it  is  continuous  on  every  compact  subinterval  of  the 
given  interval. 

8.3  Differentiation 

8.8.1  Differentiability:  If  /  ,  g  and  6  are  real-valued  functions,  and  /  and  g  are  continuous  on  a  com¬ 
pact  interval  o  ,6  with  a  <6  ,  then  we  say  g  it  a  derivative  of  f  on  o  ,6  ]  with  modulu*  of  differentiability 
6  if  and  only  if  for  all  positive  t  we  have 

!  /  (v)  -  /  (*)  -  -  *)!  ^  «  i  v  -  1  I 

whenever  z  and  y  are  in  [o,6j  and  |  y  -  z  |  <  6  (« ).  A  function  /  is  said  to  be  differentiable  on  an 
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interval  if  and  only  if  there  is  a  function  g  that  is  its  derivative  on  that  interval. 


/(»)-/  (*: 


8.3.2  Derivative:  D  =*  A (7/  )  Az  /  lim  - - - 

V  ~  x 


Comment *;  Let  g  =  Dj  f  ;  then 


/  (y )  ~  /  (* ’ 


y  -  1 

Notice  that  the  derivative  of  a  function  is  always  computable,  as  is  the  result  of  evaluating  that  derivative 
on  any  point  in  the  interval. 

8.3.3  Differentiation  Theorem:  If  /  is  differentiable  on  an  interval  7  then  D,  f  is  a  derivative  of  /  on 
that  interval.  If  /  has  two  derivatives,  then  they  are  equal. 

Proof:  Since  /  is  differentiable  on  7  we  know  that  there  g  continuous  on  7  and  there  exists  6  such 
that  for  all  e  >0  and  all  z  ,  y  in  7  for  which  |  y  -  z  |  <  S  (« )  we  have 

i  /  ( v )  -/(*)-?(*)(»-*) I  ^ 

Now,  for  any  z  in  I  let 

h(x)  =  (D,  f  )(r  )  =  /  lim  F(y), 

where  F  (y )  =  \f  (y)  -  /  (z  )]  /  (y  -  z). 

We  must  show  that  the  limit  h  (z )  exists  and  is  equal  to  g  (z  ). 

Notice  that  the  definition  of  F  includes  a  division,  which  is  computable  only  if  we  know  a  modulus  of 
nonnullity.  We  can  eliminate  this  problematic  division  be  rewriting  the  limit: 


,lim  F(y)  =*  lim  Fix  -rS[b-z) 


lim  F\x  +  (6  -z )/  n 


Now  notice  that 


F  x  (6  -z  )/  n  ]  = 


(i  -z  )/  n 


f  z  -  (b  -z )/  w  j  -  f  (x) 
b  -  z 


The  troublesome  division  by  y  -  z  has  been  converted  to  a  safe  multiplication  by  n  .  (The  division  by 
b  -  z  is  not  troublesome  because  its  modulus  of  nonnullity  does  not  have  to  vary  during  the  limiting 
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process.) 

We  now  show  that  F(y)  converges  to  g[z)  as  y  —  z.  Let  an  t  >  0  be  chosen;  we  must  find  a 
y  >  x  such  that  1  F  (y )  -  g  (x  )  |  $  t.  Take  y  <  z  i-  <5  (e ) ,  where  6  is  the  modulus  of  differentiability. 
Since  /  is  differentiable  we  know 

!/(»)-/(*)-  »(*)(v  -  *  )  I  H  t  y  -  z 

Dividing  both  sides  by  y  -  z  >0  yields 


y  -  z 


-  ?(*)!  < 


That  is,  i  F  (y  )  -  g  (z  )  |  ^  e  .  Q.E.D. 

8.8.4  Properties  of  the  Derivative:  The  derivative  of  the  sum  is  the  sum  of  the  derivatives: 

D  [sum  -  (/ )]  =  sum  °(Df,Dg) 

The  derivative  of  the  product  is  the  sum  of  the  product  of  the  First  times  the  derivative  of  the  second,  and 
the  product  of  the  second  times  the  derivative  of  the  First; 

D  [prod  0  (/  ,<7 ))  =  sum  0  [prod  °(f,Dg),  prod  •  (g  ,Df  )] 

The  derivative  of  the  reciprocal  of  a  function  is  the  quotient  of  the  negative  derivative  of  the  function 
divided  by  the  square  of  the  function. 

D  (rec ip  0  /  )  =  quo  •  [neg  •  Df  ,  (expt  2)  •  / 

The  derivative  of  the  identity  function  is  the  constant  1  function: 

D  Id  =  K  1 


The  derivative  of  a  constant  function  is  the  constant  0  function: 


D  (K  c  )  =  K  0 


8.3.5  Repeated  Differentiation:  Dn  =>  rpt  n  D 


Comments:  For  example, 


D7  f  =9  rpt  2  D  /  =*  D[D  f) 
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The  conventional  notations  ‘  dF  j  dz’  and  'DtF'  are  equivalent  to  ‘D  XzF'. 

8.S.6  Partial  Differentiation:  8  k"  =»  A  /  {A(i!  •  •  zn\[D  Xz(f  z ,  •  z4_,  z  zk  •  zn)  zt  } 

Comment i:  Thus  8*"  is  the  partial  derivative  of  an  rj-adic  function  with  respect  to  its  Jc  th  argument. 


For  example,  if  g  =82*/  then 


g(u,v,w)  =  D[\z(f  u  1  u>)j(t) 


Thus  the  conventional  notation  ‘8  F  /  8  y  where  F  is  a  function  of  z ,  y  and  z  in  that  order,  is 
equivalent  to  ‘82s  A(ryz)/”. 


8.4  Integration 


1  _  _  «-l  /  L  _  \ 

8.4- 1  Quadrature:  quad  =>  X(fabn)  -  Y)  f  la  +  i - 

l  n  i-o  l  n  ) 

Comments:  Thus  quad  /  a  b  n  is  the  quadrature  of  /  from  a  to  6  with  n  steps. 

8. 4- 2  Integral:  J  =>  A(  /  ab  )lim(quad  f  a  b)\ 

Comments :  Therefore, 

f  /  o  6  =s>  lim(quad  /  a  6  )  =s>  lim  [— - —  V  /  fa  +  t— - — 

J  n  -*ac  I  fl  4  „ n  1  H 


We  usually  write  ‘J*  /  ’  for  ‘J*  }  a  b  ’.  Thus  the  conventional  notation  F  dz  ’  is  equivalent  to 

'/V' 

5.^.5  Reversed  Range  of  Integration:  J  f  =*  /  ,  if  a  >  6  . 


8.4  4  Properties  of  Definite  Integral: 


The  integral  of  the  sum  is  the  sum  of  the  integrals: 


Ja  sum  "  (/  ,9  )  =  9um(/a  /  -  /a  ?) 


The  integral  of  a  constant  times  the  function  is  the  constant  times  the  integral  of  the  function: 


fa  prod  •  (Ka,  /  )  =  pro d(a,  f 


The  integral  from  a  to  6  is  the  integral  from  a  to  e  plus  the  integral  from  c  to  6 
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8-4-5  Indefinite  Integral: 


/:>  *  l‘>  -  ft 


Notice  that  J  f  is  an  indefinite  integral  of  /  in  the  following  sense:  If  g  =  J  /  ,  then 
9  (z  ,y  )  =  J"  /  .  Another  notion  of  indefinite  integral  is  J  f  ,  since  ifh  =  Ja/  then  h  (z)  =  fj. 

8-4-6  Properties  of  Indefinite  Integral: 

fa  sum  "(/.?)  =  sum  °  i  fa  /  -  /,?) 

Ja  prod  •  (Ka,  /  )  =  prod » (K <*<  ft  f  ) 

f.l  -  »um-  K(J//  1,  /,  /  1 

8-4-7  Fundamental  Theorem  of  Calculus: 

If  /  is  continuous  on  a  proper  interval  /  and  a  is  a  point  of  / ,  then 

D(/./>  =  / 

That  is,  the  derivative  of  the  indefinite  integral  (of  the  second  kind)  is  the  original  function.  Further¬ 
more,  if  h  is  any  differentiable  function  such  that  Dh  =  /  ,  then  for  some  real  a, 

dif»(Ja/  ,  h)  =  Ka 

That  is,  the  difference  of  J  f  and  h  is  a  constant  function. 


8-4-8  Natural  Logarithm:  In  =s»  J^recip 


Comments:  Therefore, 


In  i  =»  J^recip  *  =>  recip  =  ft  Z~'dx 


8-4-9  Definite  Integral  of  N-adic  Function: 


kjaf  =>  A(x,  -  -  I* _ i)[ Ja  Az(/  *,  •  •  •  z*_,  z  xk  i„-i)j 

Comments:  Thus  -fj  is  the  definite  integral  of  the  n-adic  function  /  ,  taking  its  k  th  argument 
from  o  to  6  .  Notice  that  the  definite  integral  of  an  n-adic  function  is  and  (n-l)-adic  function.  For 
example,  if  9  =  *  I  /  ,  then 
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c  * 

g[u,v)  =  Ja  Ax  /  (u  ,i  ,v  )j 

Thus,  the  conventional  notation  ‘J*  F  dy\  where  F  is  a  function  of  i,  y  and  z  in  that  order,  is 
r  » 

equivalent  to  J  ^{zyz  )F’. 

8-4-10  Indefinite  Integral  of  N-adie  Function: 


*■/./  -  *<*>  »*(/ 


z* -I  2  **-,1  '  '  '  M 


Comments:  As  expected,  the  indefinite  integral  of  an  n  -adic  function  is  an  n  -adic  function.  For  exam¬ 


ple,  if  g  =  2* J*  /  ,  then 


g(u,v,w)  =  f  Ax  |/  ( u  ,  x  ,w  )] 


Appendix  1:  Semlnormal  Form  of  Square-Root  of  3 

We  exhibit  the  First  few  pages  of  a  seminormal  form  for  s/2  computed  by  the  formula  derived  in  the 
body  of  the  report.  Ellipses  (  •  •  •  )  indicate  where  major  portions  of  the  formula  have  been  omitted;  how¬ 
ever,  they  look  about  the  same  as  the  portions  shown  below.  The  complete  formula  is  1974379  symbols  in 
length,  which  would  occupy  about  4700  pages  in  the  format  used  below.  It  must  be  emphasised  again 
that  the  formula,  although  large,  is  finite. 

The  formula  was  computed  by  converting  the  lambda  calculus  definitions  into  the  corresponding  LISP 
program  shown  in  Appendix  2.  The  length  of  the  formula  was  computed  in  two  ways;  (1)  by  direct 
counting  (using  a  LISP  program)  of  the  output  of  Appendix  2.  (2)  by  converting  the  Appendix  2  program 
into  the  related  length  computation  program  shown  in  Appendix  3. 

Several  recurring  features  are  visible  in  the  portion  of  the  formula  shown  here.  The  expressions 
‘(A(xy)x)’  and  ‘(A(zy)y)’  represent  the  truth  values  true  and  false.  Formulas  beginning 
‘(As  (s  •  ■  •  ))  are  list  values,  some  of  which  end  with  *j_  j_’  and  hence  represent  null  lists.  The  formu¬ 
las  ‘(A(zyz)x)’,  ‘(A(xyz)y )’  and  ‘(A(zyz)z)’  represent  the  three  basic  selector  operations  on  lists:  null,  hd 
and  tl. 


\/2  => 

(Y 

(As 

(At 

U(((Y 

(A/ 

(A(mn  ) 

((n  (A(zyz)z)) 

(A(zy )y ) 

((m(A(xyz)z)) 

(A(zy)z) 

if  (m  (A(zyz  )z  ))(n  (A(zyz  )z  )))))))) 
(A. 

(«  ( A ( zy  ) y  ) 

(As  (s  (A(zy)z) 

(As  (s  ( A ( xy  ) z  )  ^  J_  )))) 

A) 

(A(zy )z) 

((Y 

(A/ 

(A(mn  ) 

((»  (A(xyr)r)) 

(A(zy )y ) 

((m(A(zyz)z)) 

(A(zy  )z) 

(/  (m  (A(zyz  )z  ))(n  (A(zyz  )z  )))))))) 

k 

(As 

(s(A(zy)y) 

(As  (s  (A(zy  ) x  )  j_  _L  )) 

(As  (s  (A(zy)x)  j_  j_  )))))) 

(A(zy )y ) 

(A(zy  )z )) 

(As 

(«  (A(zy  )y ) 

((Ay 

(As 

(s (A(zy)y ) 
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(*  (A(xy )y ) 
y 

(A* 

(«  (A(xy)x)  X  X  ))))) 
(A* 

(x(A(xy)y) 

(A*  (*  (A(xy  )x  )  J_  J_  )) 

(A«  (x  (A(xy  )x  )  J_  —  ))))) 

(A* 

(x  (A(xy  ) y  ) 

((A» 

(A* 

(« (A(xy  )y ) 

(A(xy )y ) 

(Ax 

(«  (A(xy  )y ) 

y 

(A, 

(x  (A(zy)z) 

I  ))))))) 

( A  « 

(«(A(xy)y) 

(A*  («  (A(xy)x  )  _L  J_  )) 
(A*  (*  (A(zy )x  )  J_  X  ))))) 
(A*  (*  (A(zy)z)  X  i-  )))))) 

8 

(« (A(xy )y) 

(((A. 

(«  (A(xy  )y ) 

((((A* 

(«  (A(zy  )y  ) 

((Ay 

(A* 

(x  (A(xy )y ) 
(A(xy)y  ) 

(Ax 

(«  (A(xy)y ) 

y 

(A* 

(«(A(x 

y) 

*) 

I)))) 

(A« 

(<  (A(xy)y) 

(Ax 

(«  (A(xy)x) 

_L 

X  )) 

(Ax 

(x  (A(xy)x) 


(A« 

(«  (A(*y  )v ) 

((Ay 

(A* 

(i  ( A(xy  )y  ) 

(A(*y)y) 

(A* 

(i  (A (* 

y) 

y) 

y 

(A  i 

{•  (M* 
y) 

*) 

I  ))))))) 

(A* 

(«  (A(xy  )y ) 

(A* 

(*  (X(xy) 

*) 

X  )) 

(A* 

MM* 

y) 

y) 

(A* 

MM* 

y) 

*) 

±)) 

(A. 

(«  (A(* 

y) 

*) 

I  ))))))) 

(A* 

(*  (M*y  )* ) 

1  )))))) 

(A(*y*  )y )) 

(A(*y*  )y )) 

(A(*y  )y ) 

(A(*y  )*)) 

(A. 

(«(A(*y)y) 

((((A* 

(« (A(*y  )y ) 

((Ay 

(A« 

(*  (A(*y  )y ) 


(A(zy )V  ) 

(A* 

(«  (A(r 
y) 
y) 
y 

(A* 

(« (A(z 

y) 

* ) 

I  ))))))) 

(A* 

(t  (A(zy )y ) 

(A« 

(< (A(zy ) 

*) 

J_)) 

(A« 

(*  ( A ( x 

y) 

*) 

I  ))))) 

(A* 

(«  (A(zy)y) 

((^ 

(As 

(«  (A(z 

y ) 
y) 

(A(z 

y) 

y) 

(A* 

(i  (A(z 

y) 

y) 

y 

(A* 

(«  (A(z 

y) 

*) 

I  ))))))) 

(A, 

(«  (A(zy)y) 

(A*  (*  (A(zy  )z  )  J_  _  )) 

(A* 

(«  (A(zy  )y  ) 

(A*  (*  (A(zy  )z  )  J.  J_  )) 

(As  (#  ( A (xy  )z  )  _  _  ))))))) 

(A* 

(«  (A(zy ) 

*) 


(A(ty2)y  )) 

(A(*yz)z)) 

(A(ry2)y  )) 

(A*  (»  (A(zy  )z  )  J_  _  )))))) 

I 

(* (A(zy  )y ) 

((((Ai 

(« (A(zy  )y  ) 

((Y 

(A/ 

(A(mn ) 

((((* 

(A(zyz  ) 

*)) 

(A (zyz  ) 

y )) 

(A(zyz)z)) 

m 

((« 

(A(zyz  ) 

»)) 

(A, 

(«  (A(zy)y) 

(((((((A- 

(« (A(zy)y ) 

(((/"• 

((((«  (A{zyz)z))(A(zyz)y)) 
(A(zyz)z)) 

((Ay 

(A* 

(«  (A(zy  )y ) 

(A(zy )y  ) 

(A* 

(«  (A(*y  )y ) 

y 

(A* 

(<  (A(zy  ;Z ) 

I  ))))))) 

(A* 

(«(A(xy)y) 

(A*  (*  (A(zy  )z  )  J_  J_  )) 

(A*  (*  (A  (zy  )z)  j_  J_  ))))) 
(((n  (A(zyz  )y  ))(A(zy  )y  )(A(zy  )z  )) 

((Ay 

(A. 

(* (A(ry )y) 

(A(zy)y) 

(A* 

(« (A(ry  )y ) 

y 

(A* 

(«  (A(zy  )t) 


-  ))))))) 

((A, 

(« )y ) 

(Af  (i  (A(zy  )i)  _  _  )) 

)) 

((n  (A(xyz)z))(A(zyz)y  )))) 

((Ay 

(A* 

(« (A(xy)y  ) 

(A(zy  )x) 

(A, 

(«  (A(zy  )y ) 

y 

(A* 

(<  (A(zy  )z ) 

I  ))))))) 

(((n(A(zyz)z))(A(zyz)y )) 
(A(zyz)z )))))) 

(A(zyz)y )} 

(A(zy)y ) 

(A(zy  )z )) 

(A* 

(«  (A(zy  )y ) 

(((/"* 

((((n(A(*yz)z))(A(zyz)y)) 

(A(zyz)z)) 

((Ay 

(A. 

(»  (A(zy  )y ) 

(A(*y )y ) 

(A, 

(«  (A(zy  )y ) 

y 

(A* 

(*  (A (zy ) 

*  ) 

I  ))))))) 

(A* 

(«  (A(zy  )y ) 

(A* 

(*  (A(zy  )z) 

In 

(A* 

(« (A(zy )z) 

I  ))))) 

(((n(A(zyz)y)) 

(A(zy )y ) 

(A(zy)zj) 

((Ay 

(A* 

(* (A(zy  )y ) 

(A(zy)y ) 
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(A* 

(«  (A(zy  )V  ) 

y 

(A* 

(«  (A(z 

y) 

* ) 

I  ))))))) 

((A. 

(« (M*y)y ) 

(A* 

(<  (A(xy  )x ) 

J_ 

-L  )) 

•  )) 

((n(A(zyz)z)) 

( A(xyx  )y )))) 

((Ay 

(A* 

(i  (A(xy  )y ) 

(A(zy)z) 

(A. 

(«  (A(zy  )y) 

y 

(A* 

(* (A(z 

y) 

*) 

I  ))))))) 

(((n(A(zyr  )z)) 

(A(zyz)y )) 

(A(xyz)z )))))) 
(A(xyz)z)) 

(A(*y*  )y )) 

)))))))))) 

•) 

))))))) 

•  '  ))))) 

•  )))) 

•  '  '  ))))))) 


Appendix  J:  LISP  Program  to  Compute  Seminormal  Form  of  Square-Root  3 

(setq  false  ’(&  (x  y)  y)) 

(setq  true  ’{k  (x  y)  x)) 

(defun  Not  (x)  (list  x  false  trie)) 

(defun  Or  (x  y)  (list  x  true  y)) 

(defun  Equiv  (x  y)  (list  x  y  (Not  y))) 

(setq  lofS  ’(&  (x  y  z)  x)) 

(setq  2of3  ’(&  (x  y  z)  y)) 

(setq  3of3  '(k  (x  y  z)  z)) 

(defun  List  (n  h  r)  (list  'k  ’s  (list  ’s  n  h  r))) 

(setq  Nil  (List  true  ’undef  ’undef)) 

(defun  Cons  (h  r)  (List  false  h  r)) 

(defun  Null  (x)  (list  x  1  of 3 ) ) 

(defun  hd  (x)  (list  x  2of3)) 

(defun  tl  (x)  (list  x  3of3)) 

(defun  pair  (x  y)  (Cons  x  (Cons  y  Nil))) 

(setq  nO  Nil) 

(setq  nl  (Cons  Nil  Nil)) 

(setq  n2  (Cons  Nil  (Cons  Nil  Nil))) 

(defun  Nzero?  (n)  (Null  n) ) 

(defun  Nsucc  (n)  (Cons  Nil  n) ) 

(setq  cNsucc  (list  'k  ’n  (Cons  Nil  ’n))) 

(defun  N'pred  (x)  (tl  x)) 

(setq  rpt  (list  ’Y  (list  ’ k  ’r  (list  ’k  (f  n  a) 

(list  (Nzero?  ’n)  ’a 

(list  ’f  (list  ’r  ’f  (Npred  ’n)  ’a))))))) 

(defun  Nsum  (m)  (list  rpt  cNsucc  m)) 

(defun  Nprod  (m  n)  (list  rpt  (Nsum  m)  n  nO)) 
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(setq  cNlt  (list  ’Y  (list  ’<&:  ’f  (list  'Sc  ’(m  n)  (list  (Nzero?  ’n)  false 

(list  (Nzero?  ’m)  true 

(list  ’f  (Npred ’m)  (Npred  ’n) )))))) ) 


(defun  Nit  (m  n) 

(list  cNH  m  n)) 

(defun  Ngt  (m  n) 

(Nit  n  m)) 

(defun  Nne  (x  y) 

(Or  (Nit  x  y)  (Ngt  x  y))) 

(defun  Neq  (x  y) 

(Not  (Nne  x  y))) 

(defun  Zplus  (n) 

(pair  false  n)) 

(defun  Zminus  (n) 

(pair  true  n) ) 

(setq  pi  (Zplus  n 

1)) 

(setq  p2  (Zplus  n 

2)) 

(defun  Zmag  (n) 

(hd  (tin))) 

(defun  Zminus?  (n) 

(hd  n)) 

(defun  Zplus?  (n) 

(Not  (Zminus?  n))) 

(defun  Zzero?  (n) 

(Nzero?  (Zmag  n))) 

(defun  Zneg  (m) 

(pair  (Zplus?  m)  (Zmag  m) )) 

(defun  Zsucc  (n) 

(list  (Zzero?  n)  pi 

(list  (Zplus?  n)  (Zplus  (Nsucc  (Zmag  n) ) ) 

(Zminus  (Npred  (Zmag  n)))))) 

(defun  Zpred  (n)  (Zneg  (Zsucc  (Zneg  n)))) 

(setq  cZsum  (list  ’Y  (list  ’it  'f  (list  ’it  ’(m  n)  (list  (Zzero?  ’n)  ’m 

(list  (Zminus?  ’n)  (Zpred  (list  ’f ’m  (Zsucc  ’n))) 
(Zsucc  (list  ’f ’m  (Zpred  ’n))))))))) 

(defun  Zsum  (m  n)  (list  cZsum  m  n)j 

(defun  Zprod  (tn  n)  (list  (Equiv  (Zneg  m)  (Zneg  n)) 

(Zplus  (Nprod  (Zmag  m)  (Zmag  n))) 

(Zminus  (Nprod  (Zmag  m)  (Zmag  n))))) 

(defun  Qrat  (x  y)  (pair  x  y)) 

(setq  2/1  (Qrat  p2  pi)) 

(setq  1/2  (Qrat  pi  p2)) 

(setq  1/1  (Qrat  pi  pi)) 

(defun  Qnum  (x)  (hd  x) ) 
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(defun  Qden  (x)  (hd  (tl  x))) 

(defun  Qsum  (x  y)  (Qrat  (Zsum  (Zprod  (Qnum  x)  (Qden  y)) 

(Zprod  (Qden  x)  (Qnum  y ) )) 

(Zprod  (Qden  x)  (Qden  y)))j 

(defun  Qprod  (x  y) 

(Qrat  (Zprod  (Qnum  x)  (Qnum  y))  (Zprod  (Qden  x)  (Qden  y ) ) ) ) 

(defun  Qquo  (x  y)  (Qprod  x  (Qrecip  x))) 

(defun  Qrecip  (x)  (Qrat  (Qden  x)  (Qnum  x))) 

(defun  sqrt2  ()  (list  ’Y  (list  '&  ’s  (list  'L  'k 
(list  (Neq  nl  ’k)  1/1 

(Qprod  1/2  (Qsum  (Qquo  2/1  (list  ’s  (Npred  ’k))) 

(list ’s  (Npred  ’k))))))))) 


V 


V 
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Appendix  S:  LISP  Program  to  Compute  Length  of  Seminormal  Form  of  Square  Root  2 

(setq  false  8) 

(setq  true  8) 

(defun  Not  (x)  (add  2  x  false  true)) 

(defun  Or  (x  y)  (add  2xtruey)J 
(defun  Equiv  (x  y)  (add  2  x  y  (Not  y))) 

(setq  lof3  9) 

(setq  2of3  9) 

(setq  3of3  9) 

(defun  List  (n  h  r)  (add  2  11  (add  2  1  n  h  r))) 

(setq  Nil  (List  true  1  1)) 

(defun  Cons  (h  r)  (List  false  h  r)) 

(defun  Null  (x)  (add  2  x  lof3)) 

(defun  hd  (x)  (add  2  x  2of3)) 

(defun  tl  (x)  (add  2  x  3of3)) 

(defun  pair  (x  y)  (Cons  x  (Cons  y  Nil))) 

(setq  nO  Nil) 

(setq  nl  (Cons  Nil  Nil)) 

(setq  n2  (Cons  Nil  (Cons  Nil  Nil))) 

(defun  Nzero?  (n)  (Null  n)) 

(defun  Nsucc  (n)  (Cons  Nil  n)) 

(setq  cNsucc  (add  2  11  (Cons  Nil  1))) 

(defun  Npred  (x)  (tl  x)) 

(setq  rpt  (add  2  1  (add  2  1  1  (add  2  1  5 
(add  2  (Nzero?  1)  1 

(add  2  1  (add  2  1  1  (Npred  1)  1))))))) 

(defun  Nsum  (m)  (add  2  rpt  cNsucc  m )) 

(defun  Nprod  (m  n)  (add  2  rpt  (Nsum  m)  n  nO)) 
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(setq  c.Nlt  (add  2  1  (add  2  1  1  (add  2  14  (add  2  (Nzero?  1)  false 

(add  2  (Nzero?  1)  true 

(add  2  1  (Npred  1)  (Npred  1)))))))) 


(defun  Nit  (m  n) 

(add  2  c.Nlt  m  n)) 

(defun  Ngt  (m  n) 

(Nit  n  m)) 

(defun  Nne  (x  y) 

(Or  (Nit  x  y)  (Ngt  x  y ))) 

(defun  Neq  (x  y) 

(Not  (Nne  x  y ))) 

(defun  Zplus  (n) 

(pair  false  n)) 

(defun  Zminus  (n) 

1  (pair  true  n)) 

(setq  pi  (Zplus 

nl)) 

(setq  p2  (Zplus 

n2 ) ) 

(defun  Zmag  (n) 

(hd  (tin))) 

(defun  Zminus?  (n)  (hd  n ) ) 

(defun  Zplus?  (n) 

(Not  (Zminus?  n ) } ) 

(defun  Zzero'’  (n) 

(Nzero?  (Zmag  n) ) ) 

(defun  Zneg  (m) 

(pair  (Zplus?  m)  (Zmag  m 

(defun  Zsucc  (n) 

(add  2  (Zzero-7  n)  pi 

(add  2  (Zplus?  n)  (Zplus  (Nsucc  (Zmag  n))) 

(Zminus  (N'pred  (Zmag  n)))))) 

(defun  Zpred  (n)  (Zneg  (Zsucc  (Zneg  n)))) 

(setq  cZsum  (add  2  1  (add  2  11  (add  2  14  (add  2  (Zzero?  1)  1 

(add  2  (Zminus?  1)  (Zpred  (add  2  11  (Zsucc  1))) 
(Zsucc  (add  2  11  (Zpred  1))))))))) 

(defun  Zsum  (m  n)  (add  2  cZsum  m  n ) ) 

(defun  Zprod  (m  n)  (add  2  (Equiv  (Zneg  m)  (Zneg  n ) ) 

(Zplus  (Nprod  (Zmag  m)  (Zmag  n ) ) ) 

(Zminus  (Nprod  (Zmag  m)  (Zmag  n ) ) ) ) ) 

(defun  Qrat  (x  y)  (pair  x  y ) ) 

(setq  2  1  (Qrat  p2  pi)) 

(setq  1  2  (Qrat  pi  p2 ) ) 

(setq  1  1  (Qrat  pi  pi)) 

(defun  Qnum  (x)  (hd  x ) ) 
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(defun  Qden  (x)  (hd  (tl  x) ) ) 


(defun  Qsum  (x  y)  (Qrat  (Zsum  (Zprod  (Qnum  x)  (Qden  y)) 

(Zprod  (Qden  x)  (Qnum  y ))) 
(Zprod  (Qden  x)  (Qden  y ) )) } 


(defun  Qprod  (x  y) 

(Qrat  (Zprod  (Qnum  x)  (Qnum  y))  (Zprod  (Qden  x)  (Qden  y)))) 

(defun  Qquo  (x  y)  (Qprod  x  (Qrecip  x))) 

(defun  Qrecip  (x)  (Qrat  (Qden  x)  (Qnum  x))) 

(defun  sqrt2  ()  (add  2  1  (add  2  11  (add  2  11 
(add  2  (Neq  nl  1)  1/1 

(Qprod  1/2  (Qsum  (Qquo  2/1  (add  2  1  (Npred  1))) 
(add  2  1  (Npred  1))))))))) 
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